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The interaction between localized magnetic moments and the electrons of a one-dimensional con- 
ductor can lead to an ordered phase in which the magnetic moments and the electrons are tightly 
bound to each other. We show here that this occurs when a lattice of nuclear spins is embedded 
in a Luttinger liquid. Experimentally available examples of such a system are single wall carbon 
nanotubes grown entirely from 13 C and GaAs-based quantum wires. In these systems the hyperfine 
interaction between the nuclear spin and the conduction electron spin is very weak, yet it triggers 
a strong feedback reaction that results in an ordered phase consisting of a nuclear helimagnet that 
is inseparably bound to an electronic density wave combining charge and spin degrees of freedom. 
This effect can be interpreted as a strong renormalization of the nuclear Overhauser field and is a 
unique signature of Luttinger liquid physics. Through the feedback the order persists up into the 
millikelvin range. A particular signature is the reduction of the electric conductance by the universal 
factor 2. 



I. INTRODUCTION 

The interaction between localized magnetic moments 
and delocalized electrons contains the essential physics 
of many modern condensed matter systems. It is 
on the basis of nuclear magnetSfi heavy fermion ma- 
terials of the Kondo-lattice type^ and ferromagnetic 
semiconductors In this work we focus on the in- 
terplay between strong electron-electron interactions and 
the magnetic properties of the localized moments. Low- 
dimensional electron conductors are ideal systems to ex- 
amine this physics: The nuclear spins of the ions of the 
crystal (or suitably substituted isotopes) form a lattice of 
localized moments; these spins couple to the conduction 
electron spin through the hyperfine interaction; and the 
confinement of the electrons in a low-dimensional struc- 
ture enhances the importance of the electron-electron in- 
teractions. 

In previous work we have studied the magnetic proper- 
ties of the nuclear spins embedded in a two-dimensional 
(2D) electron gas of a GaAs heterostructure^^ and in 
13 C substituted single-wall carbon nanotubes^ (SWNTs) 
as a specific example of a one-dimensional (ID) conduc- 
tor. In this work we focus on ID more generally than in 
Ref. 0: We cover not only the case of SWNTs but also of 
GaAs-based quantum wires or different (yet not in detail 
discussed) ID conductors, under the assumption that the 
electrons are in the Luttinger liquid (LL) state as a re- 
sult of their interactions. In these systems the coupling 
between the nuclear spins and the conduction electrons 
has remarkable consequences. 

Indeed, below a cross-over temperature T* (in the mil- 
likelvin range for the considered systems) the nuclear 
spins form a spiral, a helimagnet (see Figs. [T]and[2]), 
caused by the effective Ruderman-Kittel-Kasya-Yoshida 
(RKKY) interaction induced by the electron system. 
The ordered nuclear spins create an Overhauser field that 
acts back on the electron spins. This feedback is essen- 



tial: It enhances an instability of the electron conductor 
toward a density wave order, and the electronic states are 
restructured. A gap appears in one half of the low-energy 
modes and leads to a partial electron spin polarization 
that follows the nuclear spin helix. The gap can be in- 
terpreted as a strong renormalization of the Overhauser 
field, and so as a strong renormalization of the hyper- 
fine coupling constant for the gapped collective electron 
modes. The remaining gapless electron modes in turn 
further strengthen the RKKY coupling between the nu- 
clear spins. The transition temperature T* of the nuclear 
spins can therefore lie much above the temperature that 
would be found without the feedback (called Tq below). 
In SWNTs, for instance, the feedback leads to an en- 
hancement of T* by about four orders of magnitude. 

This means that below T* there is a temperature range 
where the nuclear order and the electron order exist only 
through their mutual stabilization. The nuclear spins 
and the electrons form a combined ordered phase, even 
though the energy and time scales in both systems dif- 
fer by orders of magnitude. We remark that the or- 
der is unstable in the thermodynamic limit due to long- 
wavelength fluctuations. For any realistic sample length 
L, however, those fluctuations are cut off, the order ex- 
tends over the whole system, and T* is in fact indepen- 
dent of L. 

We discuss this physics here specifically for 13 C 
SWNTs and GaAs quantum wires because such sys- 
tems have become available for experiments recently: 
SWNTs with a purity of 13 C up to 99% have been re- 
ported in Refs. [TflUlllil 

The cleaved edge overgrowth 
metho d 14 ' 15 has made it furthermore possible to produce 
quantum wires on the edge of a GaAs heterostructure 
with LL parameters as low as 1 ^ K c = 0.4 — 0.5. For both 
systems, we predict a feedback-generated cross-over tem- 
perature T* that lies in the range of 10 - 100 mK. 

This work is related to several studies found in the lit- 
erature: NMR experiments 1 ^ on 13 C enriched SWNTs 



FIG. 1: Illustration of the helical nuclear magnetism of the single wall 13 C nanotube (SWNT), which is triggered by the RKKY 
interaction over the electron system (not shown). The nuclear spins (red arrows) order ferromagnetically on a cross-section of 
the SWNT and rotate along the SWNT axis with a period Tt/kF = Af/2 in the spin xy plane (chosen here arbitrarily orthogonal 
to the SWNT axis). The blue ribbon is a guide to the eye for the helix. The feedback of this nuclear magnetic field strongly 
renormalizes the electron system through the opening of a partial gap, due to a strongly renormalized Overhauser field, and 
so modifies the RKKY interaction. Through this strong coupling of electron and nuclear systems the combined order persists 
up into the millikelvin range. As a particular consequence the electric conductance of the SWNT is reduced by a factor of 
precisely 2. 

FIG. 2: Illustration of the helical nuclear magnetism for a GaAs-based quantum wire. Compared with the SWNT (Fig. [TJ) 
the number of ferromagnetically locked spins on a cross-section through the wire is much larger, and the wavelength it/Icf of 
the helical rotation along the wire is longer. The feedback effect remains otherwise the same and the combined electron and 
nuclear order persists into the millikelvin range as well. 



grown inside regular SWNT q 10 ' 11 have revealed the exis- 
tence of a large gap of about 30 K in the spin response. 
While the microscopic origin of this gap seems still be un- 
resolved, the NMR response could be well modeled with a 
partially gapped Tomonaga-Luttinger mode h 18 ' 19 Inter- 
estingly our microscopic theory predicts a spin excitation 
gap for a part of the transverse electronic susceptibility, 
although we consider isolated SWNTs and obtain a gap 
with a smaller magnitude below 1 K. 

The coupling between a quantum spin chain and a LL 
was studied in Ref. I20I . and it was shown that this system 
can acquire gaps as well. Such a system is very different 
from our model in that it involves a single chain of small 
quantum spins with anisotropic coupling to the electrons. 
Such an anisotropy appears spontaneously in our case, 
and built-in anisotropy has a very different effect as dis- 
cussed in Sec. IVIDI A spin gap also appears in a LL in 
the presence of Rashba spin-orbit interactions^^ 2 - LLs 
with a gap in the spin sector are known as Luther-Emery 
liquids f2£ and the partially gapped LL in our model has 
indeed a strong resemblance to such a system. Yet the 
gap opens not only in the spin sector but involves a com- 
bination of electronic spin and charge degrees of freedom, 
therefore, in addition breaks the usual spin-charge sepa- 
ration of a LL. The RKKY interaction at zero tempera- 
ture was calculated for LLs in Ref. [2H and for the case 
including Rashba spin-orbit interactions in Ref. [25l . The 
use of the hyperfinc interaction of 13 C to couple spin and 
valley quantum numbers in carbon-based quantum dots 
was explored in Ref. [26|. 

Very recent spin blockade measurements^ on quantum 



dots formed by 13 C SWNTs suggest that the hyperfine 
interaction constant Aq is by 10 2 - 10 3 larger than what 
is expected from 60 C data^l or band structure theory.— 
However, this interaction strength is inferred from the 
comparison with models that were originally designed for 
GaAs quantum dots, and so the precise value of Aq re- 
quires further investigation.— 

The observation of LL physics has been reported for 
various ID conductors such as carbon nanotubeS ) 30 ' 31 i 32 
GaAs quantum wires ^^i 3 ^ 3 ^^ bundles of NbSe 3 
nanowires^ polymer nanofibers^ atomic chains on in- 
sulating substrates^ MoSe nanowires^ fractional quan- 
tum Hall edge states^ and very recently in a bulk ma- 
terial, conjugated polymers at high carrier densities^ If 
there is a coupling to localized magnetic moments, we 
expect that the effect described in this paper should be 
detectable in these systems as well, with the exception 
of the chiral LLs of fractional quantum Hall edge states 
because they lack the backscattering between left and 
right moving modes that is crucial for the effect. To 
overcome this restriction, two edges with countcrpropa- 
gating modes would have to be brought close together by 
a constriction. 

Recently much progress has been made in 
producing and tuning the properties of car- 
bon nanotubes as quantum wires or quantum 
dots ^Q- 11 - 12 - 13 - 30 - 43 - 44 - 45 - 46 ' 47 - 48 ' 49 - 50 - 51 - 52 - 53 - 54 - 55 - 56 - 57 - 58 - 59 - 60 

and ultraclean SWNTs are now available. 5 ^ 60 

The outline of the paper is as follows: In the next 
section we state the conditions for the discussed physics 
and present the main results. A detailed account of the 
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theory is then given: In Sec. IIIII we derive the effective 
model. The nuclear order and its stability without the 
feedback is discussed in Sec. IIVI The feedback and its 
consequences are examined in Sec. [V] In Sec. I VII we 
discuss the self-consistency of the theory. The effect of 
the renormalization above the cross-over temperature is 
outlined in Sec. IVIIl In Sec. I Villi we show that the 
single-band description we have used in the preceding 
sections is appropriate for SWNTs, which normally re- 
quire a two-band model. We shortly conclude in Sec. 
IIXI The Appendices contain the technical details. The 
numerical parameters we use and derive for the SWNTs 
and GaAs quantum wires are listed in Table [IJ For a brief 
overview, we refer the reader to Sec. HU 



II. CONDITIONS AND MAIN RESULTS 

We summarize in this section the conditions and the 
main results of our work. This allows us to give an 
overview of the physics to the reader without going into 
the technical details and conceptual subtleties. These are 
then discussed in the subsequent sections. 

Two conditions for the described physics are impor- 
tant: A ID electric conductor in the LL state confined 
in a single transverse mode (in the directions perpendic- 
ular to the ID conductor axis), and a three-dimensional 
(3D) nuclear spin lattice embedded in this ID conductor. 
Higher transverse modes are split off by a large energy 
gap A t . The coupling of the nuclear spins to the elec- 
trons is weighted by the transverse mode, which even- 
tually leads to a ferromagnetic locking of the N± ^> 1 
nuclear spins on a cross-section in the transverse direc- 
tion. Consequently these fcrromagnctically locked nu- 
clear spins behave as a single effective large spin, allow- 
ing us to use a semiclassical description (see Sec. MI Bl 
and IIII C[) . This picture can be different for strongly 
anisotropic systems, where the coupling to the electron 
spin favors a different spin locking (see Sec. IVI Dp . How- 
ever, the physics described here remains valid as long as 
this locked configuration has a nonzero average magneti- 
zation. 

In addition, we treat only systems in the RKKY (Ref. 
l6l[ ) regime, which is indeed the natural limit for the 
electron-nuclear spin coupling. This regime is charac- 
terized by energetics such that the characteristic time 
scales between the slow nuclear and the fast electron dy- 
namics decouple. This makes it possible to derive an 
effective instantaneous interaction between the nuclear 
spins, the RKKY interaction, which is transmitted by 
the electron gas. If Aq is the hypcrfinc coupling constant 
between a nuclear spin and an electron state localized at 
the nuclear spin site, and if Ep denotes the typical en- 
ergy scale of the electron system, the RKKY regime is 
determined by the condition A /Ep <C 1. This condi- 
tion is naturally fulfilled in GaAs-bascd low-dimensional 
conductors^ 2 - where Aq/Ef ~ 10 -2 (for both Ga and 
As ions) and in carbon nanotubc systems grown entirely 



from the 13 C isotope (which has a nuclear spin / = 1/2) 
wher o 27 ' 28 Aq/Ef ~ 10~ 6 . Recent measurements on 13 C 
nanotube quantum dotsi 2 - suggest a much higher value 
though, but still such that Aq/Ef < 1CP 3 . An adjust- 
ment of this value, however, might be necessary because 
it relies on models that were not specifically tailored for 
13 C nanotubes^ To clarify this discrepancy between the 
reported values of Aq further experimental and theoret- 
ical work is required. We can speculate though that a 
similar renormalization of Aq as presented in this paper 
can also occur for the quantum dot system, and hence 
mimic a larger value of Aq. Due to this, the band struc- 
ture value Aq/Ef ~ 10~ 6 for the bare, unrenormalized 
hyperfine interaction strength is used in this work. 

The RKKY energy is minimized when the nuclear spins 
align in the helimagnetic order. Through the separation 
of time scales and due to the large effective nuclear spins, 
this order can then be treated as a static nuclear mag- 
netic field acting on the electrons. Most remarkably, this 
interaction is relevant in the renormalization group sense 
for the electron system and leads to the opening of a 
gap in one half of the electron excitation spectrum. This 
gap can be interpreted as a strong increase of the nuclear 
Ovcrhauscr field in the direction defined by the nuclear 
hclimagnct, while the hyperfine coupling in the orthog- 
onal directions remains unrenormalized. The gap in the 
electron system is the result of the strong binding of col- 
lective electron spin modes to the nuclear magnetization. 
The resulting RKKY interaction is then mostly carried 
by the remaining gapless electron modes and becomes 
much stronger. This leads to a further strong stabiliza- 
tion of the nuclear hclimagnct. Through this feedback 
the combined order remains stable up into experimen- 
tally accessible temperatures (see below). 

The strong renormalization is in fact possible due to 
an instability of the LL toward a density wave order^ 3 - 
which is signaled by the divergence of the electron sus- 
ceptibilities at momentum 2k p- The same divergence 
is responsible for ordering the nuclear spins, and so the 
back-action of the Overhauser field on the electrons en- 
hances the instability for a part of the electron degrees of 
freedom. This results in the partial order in the electron 
system. Due to this, the effect of the feedback is strong 
even for very weak Aq. 

We emphasize that this feedback is a pure LL effect 
and absent in Fermi liquids. It leads to a number of 
experimental signatures (described below) that may be 
used to unambiguously identify a LL without the need of 
fitting power laws to measured response functions. Let 
us also mention that an alternative test of the LL theory 
has been proposed for strongly interacting ID current 
rectifiers^MS Here a pure LL signature is found in 
form of a specific asymmetric bump in the I — V curve. 

Table Q] lists the physical parameters we use for the nu- 
merical estimates for the GaAs quantum wires and the 
13 C SWNTs. With these values we find that the feedback 
effect is most remarkable for the 13 C SWNTs where the 
cross-over temperature for the nuclear hclimagnct with- 
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Physical quantity 


GaAs Quantum Wire™ 


13 C Single Wall Nanotube 6 


Hyperfine (on-site) coupling constant |j4o| 


1 K ; 90 jtteV 


7 mK ; 0.6 fj,eV 


Nuclear spin / 


3/2 


1/2 


Electron spin S — 1/2 






Fermi vector fcj 


1 x 10 8 m" 1 


4 x 10 8 m" 1 


Fermi wavelength A_f = 2-k/Uf 


63 nm 


17 nm 


Electron density n e i = 2kp/n 


0.6 x 10 8 m _1 


2.4 x 10 8 m _1 


Fermi velocity vf 


2 x 10 5 m s _1 


8 x 10 5 m s" 1 


Fermi (kinetic) energy Ef = hvFkp/2 


7 meV 


0.1 eV 


Lattice spacing a 


5.65 A 


2.46 A 


Nuclear spin density (ID) m = 1/a 


1.8 x 10 9 m" 1 


4.1 x 10 9 m" 1 


Electron fraction per nuclear spin n e i/ni 


0.04 w 1/28 


0.06 » 1/17 


System length L 


2 - 40 ion 


2 (im 


Number of sites in transverse direction N± 


~ 50 x 50 


~ 50 


Luttinger liquid parameter (charge) K c 


0.5 


0.2 C 


Luttinger liquid parameter (spin) K s 


1 




Approximate bandwidth A a = hvp/a 


0.23 eV 


2.1 eV 


Longitudinal level spacing Al 


3-70 


260 fieV 


Transverse level spacing (subband splitting) At 


> 20 meV 


0.65 eV 


Exponent g — g XtV . z (single-band expression) 


0.75 


0.8^ 


Exponent g = g x ,y,z (SWNT 2-band expression for Tq) 




0.8^ 


Exponent g' = g' x>y 


0.67 


0.33 3 


Exponent g' z 


0.33 




Cross-over temperature (without feedback) T * 


53 mK ; 5 ^teV 


2 fiK ; 0.2 neV 


Cross-over temperature (with feedback) T* 


75 mK ; 7 fieV 


11 mK ; 1 /ieV 


Renormalized hyperfine coupling constant (in the direction 






of the nuclear spin polarization) A* = Aa^/a) 1-9 


~4.6K;~ 400 jieV 


~ 0.25 K ; ~ 22 ^eV 


Upper bound for cross-over temperature B* — SA*n e i/ni 


~ 80 mK ; ~ 7 peV 


~ 10 mK ; ~ 1 fj,eV 


Electron spin polarization (fraction) (S) / S 


0.05 


5 x 10" 5 


Electron spin polarization (on-site) {Si)/S = (n e i/nj){S)/S 


0.002 


3 x 10~ 6 


Correlation length for gapped electrons £ 


£ = £oo = 0.2 /xm 


£ = L = 2 fim 



"From Refs. [Ii[l5|[l^|35|[67i 

6 From R cfs. 27 .28.68.69 70 71. 

' Sec Sec. rvTTTl for the use of K C ,K S within the 2-band description 
of SWNTs. 



TABLE I: Physical parameters for GaAs Quantum Wires and C single wall nanotubes used in this paper, and the derived 
quantities discussed in the text. Note that the values of the derived quantities have an 0(1) uncertainty, which is unavoidable 
in any Luttinger liquid theory due to the required ultraviolet cutoff [see also the discussion before Eq. (|55jl ]. Within this 
uncertainty we have fcsT* ~ B* . Since £ = L for SWNTs, increasing L (within limits) also increases A* and B* . Anisotropic 
hyperfine interactions can slightly reduce A* and B* (see Sec. IVI D|) . Energy is converted into temperature by the relation 1 
eV = 11604.5 K. 



out the feedback, Tq [Eq. would be close to a 

microkclvin. Through the feedback, however, Tq is re- 
placed by the correct T* [Eq. ([71]) ]. which reaches into 
the experimentally accessible millikelvin temperatures^ 

For GaAs quantum wires the effect on the cross-over 
temperature is much less pronounced due to a less dra- 
matic modification of the LL parameters. Yet through 
the larger ratio Aq/Ef we already have Tq ~ 50 mK, 
which increases through the feedback to T* ~ 80 mK. 
Figures [3] through [5] show the dependences of these tem- 



peratures and of the nuclear magnetization on the varia- 
tion of different system parameters. Note that the large 
values of T * and T* for the GaAs quantum wires are due 
to the small value of K c = 0.5 we use from Ref. [l6| for 
high quality quantum wires. The more common quan- 
tum wires with weaker electron-electron interactions and 
K c ks 0.8 lead to T* - 1 mK, as shown in Fig. M In this 
figure we also show the energy scale B* [Eq. (f?4")) ]. which 
we use as an upper bound for fcsT*, below which our ap- 
proach is controlled. Further self-consistency conditions 
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arc discussed in Sec. IVII 

As a general rule, stronger electron-electron interac- 
tions (i.e. smaller LL parameters K c ) lead to larger T* 
in a much more pronounced way than a larger value of 
the hyperfme constant Aq. The explicit dependence can 
be read off from Eq. (fTTj) . and is given by 



k B T* =I\A \D' 



A a 
I\Ao\ 



(1) 



where k B is the Boltzmann constant, / the nuclear 
spin, A Q = hvp/a is on the order of the bandwidth 
(with vf the Fermi velocity and a the lattice constant), 
D' [Eq. (172")) ] a nonuniversal dimcnsionlcss constant of 
about D' ~ 0.2 — 0.3 for the values in Tabic HI and 
g' = 2K C /[K S (K C + K- 1 )] [Eq. (J68|)], where X c ,i^ s 
are the LL parameters associated with charge and spin 
fluctuations, respectively. Note that SWNTs require a 2- 
band description and so four different LL parameter o 68 ! 69 
{Kc = K cS = 0.2, K cA = K sS = K sA = 1; see Sec. 
IVIIip . While wc take this into account when we neglect 
the feedback, we show in Sec. IVIIII that the single-band 
description with K c = 0.2, K s = 1 is quantitatively valid 
when the feedback is taken into account, and therefore 
can be used for the determination of T* , and the renor- 
malized hyperfme constant A* below. 

We stress that T* is independent of the system length 
L for any realistic sample (provided that we have L 3> 
kp 1 such that the LL theory is applicable). For very large 
L there would be a cross-over where T* is replaced by a 
L-dcpcndent quantity such that T* — > as L — > oo. This 
cross-over, however, occurs only at L that lie orders of 
magnitude above realistic sample lengths (sec Sec. IIVBI 
and Appendix |P|l . 

The order parameter for the nuclear hclimagnet is the 
2kp Fourier component of the magnetization, which has 
close to T* the behavior of a generalized Bloch law [Eq. 

GOD] 



m 2 k F 



T 

T* 



3-2g' 



(2) 



This magnetization may be detectable by magnetic sen- 
sors with a spatial resolution smaller than the period of 
the helix Af/2 ~ 10-30 nm such as, for instance, mag- 
netic resonance force microscop y 72 i 73 

The nuclear spin ordering acts back on the electron 
system and leads to a strong renormalization of the hy- 
perfme interaction between the nuclear spins and a part 
of the electron modes. We can capture this renormaliza- 
tion by the replacement Aq — > A* of the hyperfme con- 
stant. We emphasize though that this replacement also 
requires a reinterpretation of the role of A* : It no longer 
describes the local coupling between a nuclear spin and 
an electron at a lattice site, but the coupling of a nuclear 
spin to a fraction of the collective electron modes in the 
LL. The modified A* has thus a similar interpretation 
as the dressing of impurity scatterin g 74 ! 75 in a LL that 



no longer corresponds to the backscattering of a single 
particle but to the generation of collective density waves 
near the impurity site. The renormalization is expressed 
by [Eq. flST" 



.4* 



(3) 



where g = (K c + K~ 1 )/2 and £ = min(L, At, Coo) is a 
correlation length. Here At = hvp/ksT is the thermal 
length, and £oo = a(A a / IAom2k F ) 1 ^ the correlation 
length for an infinite system. We stress that that £ cannot 
exceed L or At- An uncritical use of £ = £oo exceeding L 
or At can lead to self-consistency violations of the theory 
as explained in Sec. IV Al Note that for noninteracting 
electrons (including Fermi liquids) g = 1 and so Aq re- 
mains unrenormalizcd. The increase Aq — > A* is a direct 
consequence of LL physics. 

For the systems under consideration we have I\A*\ 3> 
ksT* and hence the electron spin modes following the 
nuclear hclimagnet (fcrromagnctically or antifcrromag- 
netically according to the sign of ^4o) are pinned into a 
spatially rotating spin density wave. This affects, how- 
ever, only one half of the low-energy electron degrees of 
freedom. The remaining electron spins remain in their 
conducting LL state. They do no longer couple to the or- 
dered nuclear spins, yet couple to fluctuations out of the 
ordered nuclear phase with the unrenormalized hyperfme 
coupling constant Aq, Those conducting electrons carry 
then the dominant RKKY interaction, which has a mod- 
ified form leading to the stabilization of the combined 
order up to the renormalized temperature T* . 

On the other hand, the gapped electrons have an ex- 
citation gap that is directly given by the renormalized 
Overhauser field, given by B* y = Im,2k F A* /2. There- 
fore, the nuclear magnetization m<ik w can be directly 
determined by measuring the electronic excitation gap 
through, for instance, tunneling into the systemi 15 i 16 i 35 i 36 
Since A* itself can depend on mik F through the correla- 
tion length ^oo, we have [Eq. (f59|) ] 



m 2 k f 



l 2k F 



i-(£ 



3-2 ff ' 



3-2 ff ' 



3-3 
2-3 



for m 2 k F < m 



otherwise, 



(4) 

where is such that £ = Coo for all T^fc^, > moo. 
Notice that for SWNTs, we have moo > 1 and so we 
always have B* y cx mik F - If, however, moo < 1, the 
cross-over between the two different scaling behaviors of 
B* y can be tuned by varying either L or At, depending 
on which one is smaller. 

Experiments detecting this phase can rely on two more 
effects. First, the freezing out of one half of the conduct- 
ing channels of the electron system leads to a drop of 
the electrical conductance by precisely the factor 2 (see 
Refs. d and l76h. Second, the breaking of the electron 
spin SU(2) symmetry through the spontaneous appear- 
ance of the nuclear magnetic field leads to the emergence 
of anisotropy in the electron spin susceptibility (see Ref. 
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T [K] 

FIG. 3: Magnetization m,2k F as a function of T for 13 C single 
wall nanotubes. Dashed line: without the feedback [Eq. (|38[) ]; 
solid line: including the feedback [Eq. (|70Jl]. The vertical 
lines mark the cross-over temperatures written next to them 
[Eqs. (O and iffljl]. T * is evaluated with the 2-band model 
of Sec. IVIIII while we use the quantitatively correct 1-band 
description for T* . Note that compared with Ref. ||we use a 
larger Aq (according to the results of Ref. [28h and so obtain 
a slightly larger T* . On the other hand T * is smaller because 
the exponent g is larger the 2-band model. 




FIG. 5: Cross-over temperature T* [Eq. (|71[l ] as a function 
of the hyperfine coupling constant Aq for 13 C single wall nan- 
otubes. T* follows a power-law T* cx Aq 3 3x ' v ^ — A®' 86 , and 
is plotted up to the self-consistency limit T* « vf/LI^b = 3 
K. The values about Ao ~ 10 -4 eV correspond to those de- 
duced in Ref. Il2l . Note that through the whole range we have 
k B T* ~ B* [Eq. (HJ)]. 



III. MODEL AND EFFECTIVE MODEL 




FIG. 4: Magnetization m,2k F as a function of T for GaAs 
quantum wires. Dashed line: without the feedback [Eq. (|38[) ]; 
solid line: including the feedback [Eq. (|70|l ]. The vertical 
lines mark the cross-over temperatures written next to them 
[Eqs. (|39p and (|7ip]. Since the exponent gr changes only little 
through the renormalization (g — » g'), the values of Tq and 
T* are of the same order of magnitude. 

l76l and Appendix [5} . The susceptibilities are defined by 
Eq. (fl"5]) and evaluated in Appendix From Eqs. (|A25|) 
and (|A26j) we find that for momenta q close to 2/cp and 
at T = 

X ^y{q)^{q-2k F r 2 ^'\ (5) 
Xz ( 9 )cx( (? -2A :f )~ 2 ( 1 -^), (6) 

with.^ = (K C K- 1 + K C K S )/2(K C + K- 1 ) [Eq. ©]. For 
A' s = 1 we have = <//2. At temperatures T > 0, these 
power-law singularities are broadened at q < it/ At [Eq. 
(|A16[) ] . The qualitative shape of these susceptibilities is 
shown in Fig. [7] 



A. Model 

We consider a system of conduction electrons and nu- 
clear spins expressed by the Kondo-latticc type Hamilto- 
nian 

H = HlP + A ^ ■ L < + E vtf 1 ? 1 ?- ( ? ) 

We have chosen here a tight-binding description, where 
the indices i, j run over the three-dimensional (3D) lattice 
sites Vi,Tj of the nuclear spins, with lattice constant a. 
The hyperfine coupling between the nuclear and electron 
spin on site i is expressed by the constant Ao, the electron 
spin operator S; = (Sf , Sf , Sf), and the nuclear spin op- 
erator 1, = (If, IV, If). For GaAs we have 1 = 3/2 and 
for 13 C this spin is I = 1/2. We shall generally set h = 1 
in this paper and reintroduce it only for important re- 
sults. We assume here an isotropic hyperfine interaction. 
The case of anisotropy is discussed in Sec. IVID1 

The Hamiltonian H^f describes the ID electrons (con- 
fined in a single transverse mode) and is given in detail 
below. In addition to the transverse confinement, we as- 
sume that the ID system has a length L on the order 
of micrometers that may be the natural system length 
or be imposed by gates (see Table HJ. In contrast to the 
usual Kondo-lattice model, H^P contains the here crucial 
electron-electron interactions. 

The last term in Eq. (|7|) denotes the direct dipolar 
interaction between the nuclear spins, or for the terms 
with i = j the quadrupolar splitting of the nuclear spins 
(a, P = x,y, z). Keeping this term would make the anal- 
ysis of this model cumbersome as we would have to solve 
a full 3D interacting problem. Yet those interactions 
are associated with the smallest energy scales in the sys- 
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FIG. 6: Cross-over temperatures T"o* [dashed line, Eq. ([39|l ]. 
T* [full line, Eq. iffljl]. and the bound B* [dash-dotted line, 
Eq. (|74|l ] for GaAs quantum wires as functions of the inter- 
action strength, expressed by the Luttinger liquid parameter 
K c (keeping K a — 1). The noninteracting limit is K c — 1, 
smaller K c > correspond to increasingly stronger repulsive 
interactions. The dotted line is the continuation of T* be- 
yond the energy scale set by B* , for which the validity of the 
theory becomes uncertain. The curves of T* and Tq cross 
at K c « 0.6 because the RKKY interaction J^k F defining T* 
has a prefactor that is by 1 /2 smaller than in the case without 
the feedback. Since nuclear spin order unavoidably leads to 
the feedback, T* defines the cross-over temperature for the 
order even when Tq > T* . Note that close to K c = 1, Eqs. 
(|39|) and (|71|) diverge because the cutoff S in Eq. (lAHl) was 
neglected in the further evaluation of the RKKY interaction. 
This is valid for K c smaller than ~ 0.8. Reintroducing the 
cutoff self-consistently close to K c = 1 regularizes the diver- 
gence and leads to the displayed curves (see Appendix I A Q . 




FIG. 7: Sketch of the RKKY interaction Jg [Eq J23}] or 
equivalently the spin susceptibility Xcc(q) [Eq. (|A16[) ], 



tcm. The dipolar interaction has been estimated to be 
on the order o& 10 -11 eV ~ 100 nK. For all ions consid- 
ered here, the quadrupolar splitting is absent in 13 C and 
is otherwise the largest for As with a magnitud o 77 i 78 i 79 
~ 10~ 10 cV ~ 1 /jK. These scales arc overruled by the 
much stronger effective RKKY interaction derived below 
and, in particular, are much smaller than the tempera- 
tures we consider and that are experimentally accessible. 
This allows us to neglect the dipolar and quadrupolar 
terms henceforth. 

The model ([7]) does not yet contain the confinement 
of the electrons into a ID conductor. Since we neglect 
the dipolar interaction we can focus only on those nu- 



clear spins that lie within the support of the transverse 
confining mode. This leads to a first simplification that 
is considered right below. 

B. Confinement into ID 

In this work we consider only conductors, in which 
the electrons are confined in a single transverse mode 
4>j_ . Higher transverse harmonics are split off by an en- 
ergy A t . For SWNT, this transverse level spacing is de- 
termined by§£ At = 2ttvfH/\C\, where C is the chiral 
vector describing the wrapping of the nanotube. For a 
C = {n,n) armchair nanotube we have |C| = inacc 
(with acc = a*/V3 the distance between carbon ions), 
and so A t = (27r/v3)A a /n. Due to the armchair struc- 
ture, there are 4 nuclear spins per 3acc on the cross 
section so that the number of sites on the cross-section is 
N± = An. For our choice N± sw 50 we then find n « 12, 
leading to A t as 0.65 eV. For the cleaved edge overgrown 
GaAs quantum wires, the transverse subband splitting 
has been reported-^ to exceed 20 meV. These large val- 
ues allow us to focus on the lowest transverse mode (sub- 
band) only. 

It is then advisable to switch from the 3D tight binding 
basis into a description that reflects the confinement into 
the lowest transverse mode: Let <j) ta ,<f> tl7 ■ ■ ■ label a full 
set of 2D orthonormal transverse single-electron wave- 
functions such that (j) ta = <fr±_ , and let </>j be the 3D tight 
binding Wannier functions centered at lattice sites . We 
decompose the position vector r into a longitudinal part 
r\\ and the 2D vector along the transverse direction rj_. 
Similarly we decompose the 3D lattice index i into the 
parts i|| and i±. We then write 4>i{v) = cf>i ]] (r^)(j) i± (rj_), 
and perform the basis change of the electron operators 

A.a = (to =1- \iM n ,x,c + (ti\ix)c\ utliC + ... (8) 

with a =t, |, Cj,. t n the electron operators corresponding 
to longitudinal coordinate i|| and transverse mode i„, and 

(t n \i±) = [ dr ± <j>* tn (r ± )fc x (r ± ), (9) 

for normalized wave functions <fi tn and fa^ . With the 
condition that the electrons are confined in the 4>t = <P± 
mode, averages over the operators cj tn a vanish for 
n > 1. This allows us to drop those operators from 
the beginning, and to use the projected electron oper- 
ators cj„ = (_|_ |z_i)c^ i . The 2D Wannier wave- 

1,(7 \ I - 1 - / 2||,_L,(T 

functions <f>i x have their support over a surface a 2 cen- 
tered at a lattice site, while cf>± extends over N± sites 
in the transverse direction. The normalization imposes 
that \<t>i ± (r)\ 2 oc 1/a 2 and |0±(r)| 2 cx l/N^a 2 for r 
within the support of these wavefunctions. Consequently 
(_L \i±) = Ci ± /\ / N± with Ci ± a dimcnsionlcss constant 
that is close to 1 in the support of (f>j_ and vanishes for 



those i± where <fi± — (possible phases of the Ci ± can 
be absorbed in the electron operators cj ± a ). 

The electron spin operator is quadratic in the elec- 
tron creation and annihilation operators and we obtain 
St = Sj|| ,±Cf ± /N±. The i± dependence in the Hamil- 
tonian (J7J) can then be summed out by defining the new 
composite nuclear spins 



so that the Hamiltonian becomes 



H = H 1D = Hlf 



(10) 



(11) 



This result is remarkable in that the complicated 3D 
Hamiltonian ([7]) is equivalent to the purely ID system 
(fTTj) . describing the coupling of the ID electron modes 
to a chain of effectively large spins I^. Indeed, since 
Y] i± Cf ± = Nj_ is imposed by the normalization, the 
composite spin has a length < I < IN±. As shown 
below, the maximal alignment I = INj_ is energetically 
most favorable for the RKKY interaction, such that in 
the ordered phase the composite spin can be treated as 
an effective spin of length IN±. The prefactor 1/N± to 
the hyperfine constant Aq expresses the reduction of the 
on-site hyperfine interaction by spreading out the single- 
electron modes over the N± sites. 



C. Interpretation of the composite nuclear spins I, 

It is important to stress that N± is large: For SWNT, 
N± denotes the number of lattice sites around a circu- 
lar cross-section and is on the order of N± ~ 50. GaAs 
quantum wires have a confinement of about 50 lattice 
sites in both transverse directions and so N± ~ 2500. In 
both cases N± is large enough such that the physics of 
pure (small) quantum spins does not appear. Accord- 
ingly, we shall treat the nuclear spin fluctuations below 
within a semiclassical spin-wave approach corresponding 
to an expansion in 1 / 1 Nj_ . 

We note moreover that the interaction with the elec- 
tron spin acts only on the I,., mode. Since Ci ± ~ 1 over 

most of the support of we have 1^ w J2i ± ^i\\,ix- 
Hence, all individual nuclear spins on a cross-section cou- 
ple identically to the electrons. The RKKY interaction 
acts only on I,-,, and so any nuclear order minimizing the 
RKKY interaction energy is imposed simultaneously and 
identically on all nuclear spins on a cross-section. This 
leads to the ferromagnetic locking of these nuclear spin 
shown in Fig. [5] Otherwise said, since the electrons cou- 
ple only to the transverse Fourier mode I*., describing 
the ferromagnetic alignment, any order due to the inter- 
action with the electrons can only maximize this Fourier 
component and so lead to the ferromagnetic locking. 



t 




WW 



FIG. 8: Illustration of the cross-section through the ID con- 
ductor. The N± 3> 1 nuclear spins within the support of the 
transverse confinement of the electron wave function (central 
colored region) are ferromagnetically locked and behave like 
a single large spin I = N±I. The nuclear spins outside this 
support do not interact with the electron system and are gen- 
erally disordered because their direct dipolar interaction is 
very weak. 



For SWNT, where d ± = C through rotational invari- 
ance of <j)± on the circular cross-section, the electrons 
couple only to this ferromagnetic transverse component. 
For GaAs quantum wires deviations from the ferromag- 
netic component are concentrated at the boundary of 
the confinement described by 4>±. The nuclear spins 
in this boundary layer are more weakly coupled to the 
electron spin, with an amplitude reduced by the factor 
Cf ~ l^xfojJI 2 < 1, and hence are more sensitive to 
fluctuations. 

In the following we shall interpret Ij., mainly as this 
ferromagnetic component. For GaAs quantum wires this 
cither means that we consider only very well confined 
wires where C% x = 1 over essentially the whole transverse 
cross-section, or that we restrict our attention to only 



those nuclear spins where C, 



1, and so effectively 



slightly reduce what we call the cross-section of the wire. 

The case of an anisotropic hyperfine interaction is dis- 
cussed in Sec. IVI D[ In this case spin configurations 
different from the ferromagnetic alignment are possible. 
The described physics, however, remains valid as long 
as these configurations produce a finite nuclear magnetic 
field. 



D. Effective Hamiltonians 

From the above considerations we have seen that the 
original 3D Hamiltonian is through the confinement of 
the conduction electrons equivalent to the ID Hamilto- 
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nian which we rewrite here as 



N 



(12) 



where 



'II 



now runs over the ID sites of a ID lattice 



of length L with lattice constant a. Si = Si,. j_ is now 

a ID electron spin operator, and I; = J2i ± ^ s the 
ferromagnetic component of the Nj_ nuclear spins on the 
cross-section, as discussed before. Note that even the 
reduced An /N± remains much larger than the neglected 
dipolar interaction. 

The fact that Aq <C Ep and, in particular, Ao/N± <C 
Ep, shows that the time and energy scales between the 
electron and nuclear spin systems decouple. [A more 
thorough investigation of this condition can be found in 
Sec. IVIBI ] This characterizes the RKKY regime, in 
which the nuclear spins arc coupled through an effective 
interaction carried over the electron system. Indeed, a 
change of 1^ induces a local electron spin excitation. The 
response of the electron gas to this local perturbation is 
the electron spin susceptibility Xy, describing the prop- 
agation of the effect of the local perturbation from site 
i to site j. At site j the electrons can couple then to 
the nuclear spin lj, hence inducing the effective interac- 
tion. The strict separation of time scales implies that 
this interaction can be considered as instantaneous, and 
so only the static electron susceptibility Xiji 1 ^ ~^ 0) is 
involved in the interaction. This interaction can be de- 
rived in detail, for instance, through a Schrieffer- Wolff 
transformation followed by an integration over the elec- 
tron degrees of freedom as in Refs. OH. The result is the 
effective Hamiltonian for the nuclear spins 



H 



eff 



ij,af3 



ja/3 

J ij fa f/3 

Nl 1 



(13) 



with a, /3 = x,y, z and the RKKY interaction 



J 



af3 



^0 2 a/3 
r> " Ay ) 



(14) 



where xtf * s the static electron spin susceptibility 



a/3 
Xij 



a -10 



dte-*{[S?(t),S$(0)]), (15) 



for an infinitesimal 7/ > and the average determined 
by H e i. Note the 1/a 2 in this definition, which allows 
us to pass to the continuum limit xtf ~~* X a ^( r ) with- 
out further complication (see Appendix IA 1|) . We assume 
that the total spin is conserved in the system, and so 
X a ^{ r ) = ^a/3Xa(f)i J a P{r) = S a pJ a (r). In momentum 
space we obtain 



rreff _ J_ \ " " q fa fa 

n n TV TV 2 q q 

q , a -L 



(16) 



with TV = L/a, the Fourier transform I" = J^i e lqTi I", 



and 



A 2 A 2 r 



-tqr 



Xa[r) 



(17) 



In this derivation of the RKKY interaction we have tac- 
itly assumed that the electrons are unpolarizcd. This will 
be no longer the case once the feedback coupling between 
electrons and nuclear spins has been taken into account. 
The necessary modification to Eq. (Til)]) is discussed in 
Sec. IV El and leads to the Hamiltonian (1661) . 

This feedback is driven by the Overhauser field gener- 
ated by the nuclear spins acting back on the electrons. 
To model this we rely again on the separation of time 
scales, which allows us to treat the Overhauser field as 
a static external field for the electrons. Hence, a mean 
field description of the nuclear Overhauser field is very 
accurate. This leads to the effective Hamiltonian for the 
electron system 



rreff 
H el 



ttID 



(18) 



with hi = £ 



Ij), and where the expectation value is 
taken with respect to H^ s . The Hamiltonians H^f and 
.ff„ , and so the properties of the electron and nuclear 
subsystems self-consistcntly depend on each other. 



E. Electron Hamiltonian 

The confinement of the conduction electrons in the sin- 
gle transverse mode <f>j_ makes the electron Hamiltonian 
strictly ID. Through this dimensional reduction electron- 
electron interactions have a much stronger effect than in 
higher dimensions. In particular, they lead to a departure 
from the Fermi liquid paradigm, and often the Luttingcr 
liquid (LL) concept, based on the Tomonaga-Luttingcr 
model, is the valid starting point to characterize the sys- 
tem properties. 

We therefore consider a ID system of length L with 
electron-electron interactions that are effectively short- 
ranged due to screening by gates. Such a system allows a 
description by the Tomonaga-Luttinger model, given by 
the Hamiltonian^ 

Hi? 



= ^2 j drv F ^\ a {r)iV^ La {r)-^ Ra {r)iV^ RrT {r) 
+ Y, I drdr'U{r-r , )i>l{r)^l{r , )i> al {r , )^r), (19) 



where i\)i a (r) are the operators for left moving {1 = L = 
— , with momenta close to — fcf) and right moving (£ = 
R = +, with momenta close to +kp) electrons with spin 
(T =t= + and a =[= — . The positions r,r' run over a 
system of length L such that L ^> n/kp, where 2kp /ir = 
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n e i is the electron density of the ID conductor. The 
operator ip a ( r ) = ^La( r ) + V'-Ro-M is the conventional 
electron operator. The potential U(r — r') describes the 
screened electron-electron interaction. 

For simplicity, we consider here a single-band descrip- 
tion of the ID conductor. This is not correct for carbon 
nanotubes, which require a 2-band model. In Sec. IVIII1 
however, we show that the main conclusions and results 
are quantitatively determined by the single-band model, 
so that we can avoid the more complicated 2-band de- 
scription henceforth. 

In the Hamiltonian (fl"9)) we have assumed a linear 
fermionic dispersion relation for the left and right moving 
electrons, e q = ±v F {q T k F ) (setting the chemical poten- 
tial fj, = 0). If this assumption is valid, we have with 
the bosonization techniqu o 63 i 81 a powerful tool to evalu- 
ate the properties of the electron system to, in principle, 
arbitrary strength of the electron-electron interactions, 
leading to the LL description. 

We assume here that this theory holds. Possible devi- 
ations are discussed in Sec. IVI Al The derivation of the 
bosonic theory can then be done in the standard way^ 3 - 
by expressing the fcrmion operators ^>n a in terms of boson 
fields 6trr as 



-0£<r (r) 



(20) 



F. RKKY interaction 

The calculation of the RKKY interaction, i.e. the cal- 
culation of the electron spin susceptibility, is standard in 
the LL theory. We have outlined its derivation in Ap- 
pendices and [B] (the real space form of the RKKY 
interaction at T = has been derived before in Ref. l24h . 
and from Eq. (|Al6j) together with Eq. (fT7|) we obtain, 
for q > 0, 



Jq{g a ,v F ) 



A 2 



2-2g a 



k B T, 

r(f-^( 9 -2fc F )) 



r( 



2-5„ 



-i£(q-2k F )) 



(23) 



In this expression we have neglected an additional small 
term depending on q + 2k F and the small forward- 
scattering contribution. We note, however, that J°: q = 
Jq. Much of depends on the quantities 

g = g x , y = [K c + l/JQ/2, g z = [K c + K s )/2. (24) 

For SWNTs these definitions must be modified due to the 
existence of two bands, which we use for the case when 
the feedback to the electron system is neglected. From 
the discussion in Sec. IVIIII we have 



with r]£cr the Klein factor removing a (£, a) particle from 
the system and 



1 
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[4> c -l6 c + a{cb s -l6 s )}. 



(21) 



Here </> CjS are boson fields such that — \7(b cs \/2/n mea- 
sure the charge and spin fluctuations in the system, re- 
spectively. The boson fields 8 C S are such that V# CiS /7T 
are canonically conjugate to (b c , s . The Hamiltonian (Til))) 
can then be rewritten in these boson fields as 



ttID _ 

n„ } — 



v l ,K u (ye u {r)f 



(22) 

where K c ^ s are the LL parameters for the charge and 
spin density fluctuations, and v c s are charge and spin 
density wave velocities. The electron-electron interac- 
tions are included in this Hamiltonian through a renor- 
malization of K c s and v c iS . The noninteracting case is 
described by K c = K s = 1 and v c = v s = vf- Repul- 
sive electron-electron interactions lead to < K c < 1. 
If the spin SU(2) symmetry is preserved K s = 1 other- 
wise K s > 1. The case K s < 1 would open a gap in 
the spin secto r 63 i 81 and is not considered here. For ideal 
LLs one has w C;S = vf/K c ^ s . With Eq. (|2"2"|) we have 
furthermore assumed that kp is not commensurate with 
the lattice spacing. Altogether, this allowed us to drop 
irrelevant backscattcring and Umklapp scattering terms 
in Eq. 



9 = g x , y .z = (Kc + 3)/4. 



(25) 



We stress, however, that the single-band values (|24|) are 
quantitatively correct also for SWNTs when we take the 
feedback into account (see below). 

For K s = 1 it follows from K c < 1 that g a < 1. The 
prefactor in Eq. (|23|) is given by 



C(g) 



sin(7r 5 ) . 



T^l-g) (2tt) 



2 9 -4 



we have introduced the thermal length 

Ay = hvp/ksT, 

and the energy scale 

A a = hv F /a, 



(26) 



(27) 



(28) 



which is on the order of the bandwidth. 

A sketch of J q is shown in Fig. [7] This interaction 
has a pronounced minimum at q = ±2kp with a width 
~ 7t/At, and a depth of 



2k, 



A a 



°-C(g a ) 



\k B T) 



2-2g a 



T(ffa/2) 



r(i-W2) 



(29) 



For \q — 2k p\ > tt/Xt the zero temperature form of this 
curve becomes quickly accurate: For q > and combin- 
ing Eqs. (|A15|) and ([T7J) we find 



Al sin(7rg a ) 
A a 8tt 2 



a(q — 2k F ) 



2-2g a 



(30) 
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which can also be verified by letting T — > in Eq. ((23)) . 

The singular behavior at ±2kp defines the real space 
form J a (r). The latter can be obtained by Fourier trans- 
forming Eq. (J23J) or by time integrating Eq. (|A13j) . The 
latter is done explicitly in Appendix [Bl From Eq. (|B5P 
we then find 

_ A 2 X T sm(Trg a ) f Tza/\ T \ 2g ° 

^ W _ Stt^^ co& (^) ^ Sin h(7r|r|/A T ) J 

x F (1/2, 5q ; 1; - sinh- 2 (7r|r|/A T )) , (31) 

where F is the Gaussian hypergeometric function, de- 
fined in Eq. (|B4|) . For \r\ <C At, which is the case 
most of the time for the systems considered here, and 
for g > 1/2 we obtain the asymptotic behavior [see Eq. 
(IB71) and Ref. H 

J a (r) ~ cos(2fc F r)(a/|r|) 2 ^- 1 . (32) 

Stronger electron-electron interactions lead to smaller g a 
and so to an RKKY interaction that extends to longer 
distances. Since the order discussed below is due to the 
long-range part of the RKKY interaction, this leads to 
a better stabilization of the order and consequently to a 
higher cross-over temperature. 

Let us note that Eq. (|32|) cannot be extended to 
g a < 1/2, where it would describe an unphysical growth 
of J(r) with distance. This regime is actually not reached 
for conventional LLs with K c > and K s = 1, yet with 
the feedback below we obtain renormalized g' a < 1/2. 
Eq. (|32p then is regularized by further temperature- 
dependent corrections coming from the expansion of the 
hypergeometric function (see Appendix [Bj or at low tem- 
peratures by cutoffs such as the system length. Since the 
relevant temperatures determined below are such that 
L < At or at most L ~ At, the values g' a < 1/2 then 
lead to a RKKY interaction that decays only little over 
the whole system range. 

IV. NUCLEAR ORDER 

A. Helical magnetization 

We have seen above that the ferromagnetic locking of 
the N± nuclear spins in the direction across the ID con- 
ductor leads to a ID nuclear spin chain of composite nu- 
clear spins L with maximal size IN± S> 1. This allows 
us to treat the nuclear subsystem scmiclassically>i£ Pure 
quantum effects such as, for instance for antiferromag- 
netic chains, the Haldane gap for integer quantum spin 
chain a 82 ' 83 or Kondo lattice physics 2 - are absent, since 
their effect vanishes exponentially with increasing spin 
length. 

In the present case the starting point is the classical 
ground state of the nuclear spins described by the Hamil- 
tonian (fTS")) . The RKKY interaction J q reaches its mini- 
mum at q = ±2kp, and the ground state energy is mini- 
mized by fully polarized Li = INj_ describing a helix with 



periodicity wave vector 2kp- The corresponding ground 
states fall into two classes of different helicity, 

I' ; = IN±[e x cos(2fcTfi) ± &y sin(2/cT^i)]; (33) 

where e Ii9 are orthogonal unit vectors defining the spin 
x, y directions. Through the spontaneous selection of the 
directions e x ,y any rotational symmetry of the Hamilto- 
nian (fT3"|) is broken. 

The simultaneous existence of both hclicities cannot 
occur for these classical ground states because their su- 
perposition would result in a single helix, but with re- 
duced amplitude. A coherent quantum superposition of 
such states, on the other hand, can be excluded because 
each state involves the ordering of the large effective spins 
/ = IN± over the whole system length L, and hence 10 5 - 
10 8 individual spins. This symmetry breaking and so the 
selection of a single helicity is in fact crucial for the feed- 
back effect described in Sec. |Vj A single helix leads to a 
partially gapped electron system, while a superposition 
of the two helicities would result in an entirely gapped 
electron system. The physics of the latter case is inter- 
esting on its own, but does not occur in the present case. 

Transitions between both helical classes would involve 
a reorientation of the entire nuclear system (and through 
the feedback of the electron system as well) . Such a tran- 
sition, as well as the spontaneous emergence of domain 
walls, can be excluded because its energy cost scales with 
the system size due to the long-range RKKY interactions. 
The low-energy fluctuations about the ground state are 
either rotations of the entire nuclear spin system as a 
whole or magnons. 

Rotations of the whole system do not reduce the local 
magnetization, yet they may lead to a zero time average. 
Due to the the aforementioned separation of energy scales 
between the nuclear and electron system, however, the 
momentary nuclear spin configurations acts like a static, 
nonzero field on the electrons. Our analysis, therefore, 
is not influenced by these modes. Moreover a pinning of 
those modes, since they involve the rigid rotation of the 
entire system, is very likely. 

More important are magnons, which describe the low 
energy fluctuations to order 1/IN±. The magnon spec- 
trum for the nuclear helimagnet is derived in Appendix [Cl 
(see also Ref. HI). For the isotropic or anisotropic (when 
the feedback on the electrons is considered) RKKY inter- 
action there exists a gaplcss magnon band with dispersion 
given by Eq. (|C6[) , 

u q = 2{IN 1 _){r 2kF _ q - r 2kF )/Nl. (34) 

Let mi = (Ij) • V i /(IN±) 2 measure the component of 
the average magnetization along 1^, normalized to — 1 < 
to.; < 1. The Fourier component m,2k F then acts as an 
order parameter for the nuclear helical order. We can 
choose < TO2fe F < 1 by rotating the axes e Ii9 if neces- 
sary. Magnons reduce this magnetization as follows^ 

m2 * F = 1 -/ivTiv"?e^-l' (35) 
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where the sum represents the average magnon occupation 
number, and the momenta q run over the first Brillouin 
zone q e [— tv /a, Tv/a]. 



B. Absence of order in infinite-size systems 

In the thermodynamic limit L — > oo the sum in Eq. 
35|) turns into an integral that diverges as L/a = N 
due to the q — > magnon occupation numbers, showing 
that long-wavelength modes destabilize the long-range 
order. It is noteworthy that this is not a consequence 
of the Mermin- Wagner theorem^ The Mermin- Wagner 
theorem forbids long-range order in isotropic Heisenberg 
systems in low dimensions with sufficiently short-ranged 
interactions. An extension to the long-ranged RKKY in- 
teractions has been recently conjectured^ for the case of 
a free electron gas. The theorem thus cannot be applied 
for systems where the long-range RKKY interaction is 
modified by electron-electron interactions. Indeed, we 
have shown in previous work>£>£ that in this case long- 
range order of nuclear spins embedded in 2D conductors 
becomes possible. In the present ID case, however, the 
divergence of the magnon occupation number at q — > 
provides a direct example where long-range nuclear mag- 
netic order is impossible in the L — > oo limit. 

Realistically we always deal with samples of a finite 
length L though. The singularity at q — > is cut off 
at momentum n/L, and the q = is absent in samples 
that are not rings. This means that the sum in Eq. (|35p 
is finite, and a finite magnetization is possible at low 
enough temperatures. We shall actually see below that 
even though the cutoff at tt/L plays a significant role for 
the stability of the order in realistic systems, the mag- 
netization TO2fc F is fully determined by L-indcpcndent 
quantities. This is much in contrast to what we would 
anticipate from the L — ► oo limit. 



C. Order in finite-size systems 



o,/ L /k B T w 2I\ J$ kp \/N ± k B T = 1, we have for T > T M0 



m 2 k F (J 1 ) = 1 



1/IN A 



1 



T 



IN ± \T M0 



3-2g 



(37) 

with g = g x>y = (K c + K~ 1 )/2 from Eq. ([24]) for GaAs 
quantum wires and g = (K c + 3)/4 for SWNTs (sec Sec. 
IVIIip . We eventually write the magnetization as 



3-2g 



(38) 



where we have defined the temperature T * by 
k B T* = {IN x )^k B T M0 = 2I 2 \J* kp \ 

2 ' 2CJ r 2 ( 5 /2) 



2I 2 C(g) 



I\A \D 



A\ ( 






k B T* 


A Q 


l — 2g 




^ 3-2g 



r 2 (l-g/2) 



(39) 



where C(g) is defined in Eq. (|26p . A a = hvp/a [Eq. 
([28]) ] . and where D is the dimensionless constant 



D 



sim% 5 )r 2 (l- 3 )(27r) 



2g-i_ 



r 2 ( ff /2) 



r 2 (i- 5 /2) 



(40) 

Eq. (|38|) can be considered as a generalized Bloch law for 
the nuclear magnetization with an exponent 3 — 2g that 
depends on the electron-electron interactions. 

The arguments above are based on the assumption 
k B T < Al so that Eq. ([38]) is in principle only valid 
if k B TX <C A^. In Appendix [Dl we show, however, that 



Eq. (|38[) remains valid far into the range k B T > A^. 
In particular, it remains valid for GaAs quantum wires, 
where we find that k B T^ , k B T* ~ A^, where T* is the re- 
sulting cross-over temperature when taking into account 
the feedback onto the electron system [Eq. ([71))]. 

This means that the system length L has no influence 
on the nuclear magnetization for any realistic SWNT and 
GaAs quantum wire system. 



The energy representation of the momentum cutoff at 
tt/L is the level spacing 

A L = hv F /L. (36) 

This level spacing must be carefully compared with any 
other temperature scale characterizing the magnetization 
m 2kF of Eq. |35[) . 

In particular, if k B T < A^, the momentum quantiza- 
tion tt/L is larger than the inverse thermal length, tt/Xt- 
Since the width of the minimum of J 2 k k +q ^ s on the or- 
der of tt/Xt (see Fig. [7J, the first possible magnon energy 
u>k/l is already very large and close to the maximal value 
2/ 1 J 2 k |/iVj_. If we define a temperature Tmo at which 



D. To sets the only available energy scale 

It may seem surprising that T * does no longer de- 
pend on ./Vi. Yet we need to recall that the reduction 
of the RKKY interaction J by 1/N]_ is compensated 
through the coupling of two composite spins containing 
each N± spins /. The coupling energy therefore depends 
on {J/Nj) x {IN ±) x {IN a.) = I 2 J, which no longer con- 
tains iV_i_ (see Fig. [9|). The cross-over temperature (since 
the L dependence has been ruled out) can then only de- 
pend on the energy scales that characterize I 2 J. Since J 
is described by the width and the depth of its minimum, 
which depend both on T, there is only one characteris- 
tic temperature that can be self-consistently identified, 
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N 





























Jij/Nl 





HIP 



H 0v with HIP given by Eq. (|22[) and 



TV i 



J 



FIG. 9: Illustration of the RKKY coupling between two 
large spins composed of N± individual spins at sites i and 
j. The hyperfine interaction is reduced by distributing a sin- 
gle electron over N± nuclear spins, A = Ao/Nx, resulting 
in an RKKY interaction Jij/Nj_. This reduction by 1/Nj_ is 
compensated because N± nuclear spins are mutually coupled 
through the same RKKY interaction Jij. 



by setting I 2 \J 2 k F {T)\ = k B T. The result is T * [up to 
the factor 2 in Eq. (|39|) ]. which consequently must set 
the scale for the cross-over temperature, independently 
of the chosen approach, be it magnons (as here), mean 
field or more refined theories. 

In the next section we see, however, that T * is strongly 
renormalized by a feedback coupling between the electron 
and nuclear spin systems, which modifies the shape of 
J q itself. The feedback in addition introduces a second 
scale through a partial electron spin polarization that 
acts like a spatially inhomogeneous Zceman interaction 
on the nuclear spins. 



V. FEEDBACK EFFECTS 

A. Feedback on electrons 

We have seen that the electrons enforce a helical or- 
dering of the nuclear spins, and we have assumed that 
this helical ordering defines the spin (x, y) plane [see Eq. 
(|3"3"|) ]. In the following, we analyze the feedback of the 
nuclear spin ordering on the electrons using the Rcnor- 
malization Group (RG) approach. Since the dynamics 
of the nuclear spins is much slower than that of the elec- 
trons, we can safely assume that the main effect of the nu- 
clear spins is well captured by a spatially rotating static 
magnetic field of the form: B(r) = B x cos(2kpr)e x + 
By sin(2fci?r)e. y , with B x = B y = IAoni2k F - Note 
that we explicitly (and arbitrarily) choose the coun- 
terclockwise helicity for the helical ordering. The ef- 
fective Hamiltonian for the electron system then reads 



Hov = j drB(r)-S(r), 



(41) 



is the coupling to the nuclear Overhauser field. 

Using the standard bosonization formulas^ Hov is 
expressed as 



Hov — I - — B xy 



dr 
2ixa 



cos(V2((f> c + 9 S )) 
+ cos(V2(0 c - 6 S ) - Ak F i 



where 



B xy = IA m 2 k F /2 1 



(42) 



(43) 



and where we have not written the forward scatter- 
ing part because it has no influence. The last term is 
oscillating and is generally incommensurate except for 
Akpa = 27m, with n integer. This special case would cor- 
respond to a fine-tuning of Ep to about 1.7 eV for carbon 
nanotubes or 0.2 eV for GaAs quantum wires, which is 
quite unrealistic for the systems we consider here. We 
therefore assume in the following that Akpa ^ 2%n, and 
hence drop the last incommensurate term in Eq. (|42[) . 
The remaining cos(\/2(0c + S )) term has the scaling di- 
mension (K c + K~ r )/2. For the systems under consider- 
ation this operator is always relevant, so that the interac- 
tion term is always driven to the strong coupling regime. 
To see this more clearly we change the basis of boson 
fields by introducing the boson fields (f> + ,<j)-, 9+ , 9- , de- 
fined by22 



(44) 
(45) 
(46) 
(47) 



VkTk 

where we have set 

K = K, 




KZ 



(48) 



The new fields obey the standard commutation relations 
[0 K (x), 9 K >{y)\ = i7r<5 KjK 'Sign(a; - y) with k, k' = +, -. In 
this basis, the electron Hamiltonian reads 



H el = 



+ —^cos(V2Kfa 
a 



+ (V c - V s ) 




(49) 



(V0 + )(V0_)-(V0+)(V0_ 
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with 



v+ = (v c K c + v s K~ l )/K 
v- = {v c K~ l +v s K c )/K. 



(50) 
(51) 



If v c = v s , the electron Hamiltonian separates into two 
independent parts H e i = + HZ where is the 
standard sine-Gordon Hamiltonian, while H~, is a free 
bosonic Hamiltonian. The cosine term is relevant and 
generates a gap in the '+' sector. If v c ^- v s the terms 
in V#+V#_ and V0+V0_ are marginal and are much 
less important than the strongly relevant cosine term. In 
a first approximation, we neglect these terms. We will 
come back to this point in Sec. IV CI below. The RG 
equation for B xy then reads^ 



dl 



1 



Kr 



2 2K S 



B , 



(2 - g)B xy , (52) 



where I is the running infrared cut-off and 
g = K/2 = (K c + K~ l )/2. 



(53) 



We use this g for both GaAs quantum wires and SWNTs, 
in contrast to the g of Eq. (|2"5"|) that must be used for 
SWNTs in the absence of the feedback. As explained in 
Sec. IVIIII this is due to the fact that the feedback acts on 
each of the two Dirac cones of the the SWNT dispersion 
relation separately, and hence effectively splits the two 
bands of the SWNT into separate single-band models 
within each cone. 

Under the RG flow, B xy (l) grows exponentially as does 
the associated correlation length £ = a(l) = ac l . The 
flow stops when either £ exceeds L or At, or when the 
dimensionless coupling constant^ y(l) = B xy (l)/A a (l), 
with A a (l) = hvp/a(l) becomes of order 1. From the 
latter condition we obtain a correlation length 



£ = Coo = a[y(0)} 



-l/(2- 3 ) _ 



a[B xy /A a ] 



-1/(2-9) 



(54) 



We emphasize that with the cutoff criterion y(l) ~ 1 the 
magnitude of the resulting B xy has an O(l) uncertainty. 
In fact, we use here a different cutoff criterion as in Ref. 
d, namely A a instead of Ep. While the use of both cutoffs 
is generally justified for the perturbative RG scheme used 
here (the cutoff must be on the order of the bandwidth) , 
we notice that when using Ep we obtain for the GaAs 
quantum wires too large values for B xy that exceed Ep. 
This is unphysical as it would imply that more electrons 
are polarized than are contained in the system, and so 
it just means that B xy becomes comparable to Ep. The 
consistency of the RG scheme then requires that the RG 
flow must be stopped earlier, and the natural scale in the 
kinetic part of the Hamiltonian (|49|) is set by A a £^ Due 
to the resulting smaller gap B xy , however, for SWNTs 
the correlation length £ would exceed the system length 
L at the new cutoff scale (while £ « L in Ref. @). Hence, 
for SWNTs the flow is cut off even earlier at L. 



B. Renormalized Overhauser field and gap for 
electron excitations 



Independently of the precise form of the correlation 
length £ we can always write y(l) = y(0)(£/a) 2 ~ 9 . Since 
furthermore A a (0) / A a (l) = a/£, we obtain the following 
result for the gap, i.e. the renormalized Overhauser field 
B 



■nj 



B xy (l), 



K v = B xy (C/a) 1 " 9 , 



with 



£ = min{L, A T ,£oo = a^/AJ" 1 /^-^ } 



(55) 



(56) 



Since B xy cx Aq we can translate this directly into a 
renormalized hypcrfinc interaction constant 



A* = A (e/a) 1 ^. 



(57) 



It is important to notice that even though A* can be 
called a "renormalized hyperfine interaction" it no longer 
can be interpreted in the same way as Aq. It does not de- 
scribe the on-site interaction between a nuclear spin and 
an electron spin, but results from the reaction of the en- 
tire electron system to the ordered nuclear spin system. 
We can see this as analogous to the strong growth of 
an impurity backscattering potential in a LL ) 74 ' 75 which 
then no longer corresponds to the coupling between the 
impurity and an electron but involves a collective screen- 
ing response by the electron system. 

The values of A* are listed in Table [fl Since I A* > 
ksT, for all temperatures within the ordered phase, we 
find that precisely one half of the degrees of freedom, the 
</>_!_ fields, are gapped, while the fields remain in the 
gapless LL state. As we have shown in Ref. [7(1 this has 
the direct consequence that the electrical conductance 
through the ID system drops by the factor of precisely 2. 
Since the gap is identical to the nuclear Overhauser field, 
and so proportional to the nuclear magnetization rri2k F , 
it in addition allows to directly measure the nuclear mag- 
netization through a purely electronic quantity, the gap 
B xy . Using Eq. for m 2 k F , which shall be proved 
explicitly in Eq. (|7D|) below, we can rewrite the gap as 



IA* 



1-1^ 

T* 



3-2g' 



(58) 



with T* given by Eq. ([T]), or Eq. (n} below. For 
SWNTs A* is independent of m^kp > and so the gap B* y 
is directly proportional to the magnetization m,2k F ■ For 
the GaAs quantum wires, A* is independent of m2k F 
for small magnetizations such that £oo > L, and the 
correlation length is set by L. For a magnetization 
tnik F > n^oo = (A a /IA )(L/a) 2 ~ 9 we have £oo < L and 
A* becomes a function of m^kp- The gap therefore fol- 
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lows the curve 



Ky « 



m 2 k F = 



3-9 
, 2 ~9 



l-(^) 



,3-2g' 



1 - 



' T 
\ T* 



.3-2g' 



for m 2 fc F < m c 



otherwise. 



(59) 

Notice that the value Woo for the cross-over magnetiza- 
tion can be tuned by the system length L. If Xp < L, 
the L in the definition of moo is replaced by At- 

The physical meaning of the gapped field cj>+ is best 
seen by rewriting it in terms of the original boson fields 
4>ia using Eq. (21), 

4> + = (4> c + s )/Vk=(0 ri +^)/V2K. (60) 

A gap in the '+' sector means that a linear combination 
of the spin I electron right movers and spin j electron 
left movers is gapped. This combination is pinned by the 
nuclear helical state. This can be seen as the analog of a 
spin/charge density wave order except that it involves a 
mixture of charge and spin degrees of freedom. As shown 
in Sec. IV D[ this pinned density wave corresponds to an 
electron spin polarization following the nuclear helical 
order. 



the helix of the nuclear Overhauser field and is parallel 
(for a ferromagnetic Aq) or antiparallel (for an antiferro- 
magnetic Aq) to the nuclear spin polarization. 

Within the LL theory, the form of this polarization 
can be found very easily: The forward scattering contri- 
bution of S z oc V<^ s has a zero average. There remain 
the backscattering parts given by the Og DW operators 
in Eqs. (|A5[) - (|A7[) . Since averages over exponents 
consisting of single boson fields vanish in the LL theory, 
(e 1 *-) = (e 10 -) = (up to finite size corrections on the 
order of a/L), only the contribution to the spin density 



average consisting uniquely of the 



3 ±iV2if0a 



operators 



are nonzero because (f> + is pinned at the minimum of the 



cosine term in Eq. (49]). For Aq > this minimum is at 
V2K(j) + = -7r and for A Q < at \/2K(j) + = 0. The two 
backscattering parts depending only on <f> + arc 



(^ItWV'fllW) 



Jlikpr 



a iV2K<f> + 



2iza 



2ttq 



-sign(A)), 



(63) 

and the conjugate (ip R i (r)ipL^ (?")}• This leads to the LL 
result for the electron spin polarization density (with S = 
1/2) 



C. Corrections by the marginal terms 

The results of the RG analysis above remain almost un- 
changed if we take into account the terms V# + V#„ and 
V(/>+V(/>_. To support this assertion, we have checked 
numerically that that y = B xy /A a reaches its cutoff 
scale, while the other coupling constants remain almost 
unchanged. Following Ref. l63l . we then expand the co- 
sine term up to second order. This provides a mass term 
oc B xy for the mode. Within this approximation, the 
<j) + , 6+ bosonic fields can be exactly integrated out in the 
quadratic action. This leaves us with an effective Hamil- 

6- , which has precisely the 



tonian H cS for the fields 



same form as H el up to some irrelevant terms, 



H cil = V 



cir 



dr 
2^ 



1 



(61) 

but with the renormalizcd parameters = y/1 — T and 
v cE = v-VT^T, with 



4K S ( Vc K c + v 8 Ks"){v c K7' + v s K c ) 



(62) 



For the systems under consideration the factor \/l — T is 
slightly less than 1 and this renormalization has indeed 
no quantitative consequences. 



D. Electron spin polarization 

The pinning of the + modes leads to a partial polar- 
ization of the electron spins. This polarization follows 



(S(r)) 



LL 



-s 



sign(Ao) 



cos(2fci?r) 
sin(2fci?r) 




(64) 



The correct prcfactor, the polarization density, cannot 
be obtained from the LL theory, which only provides the 
dimensionally correct prefactor l/ira. This unphysical re- 
sult is a direct consequence of neglecting bandwidth and 
band curvature effects in the LL theory, which can lead 
to a violation of basic conservation laws. Concretely we 
obtain here an electron polarization oc 1/a which largely 
exceeds the electron density oc hp in the system. To 
cure this defect, we use the following heuristic argument 
in the spirit of Frohlich and Nabarro>i The process of 
opening the gap B* y in the </>+ field is carried mostly 
by the electrons within the interval B* y about the Fermi 
energy (using a free-electron interpretation). Hence, the 
polarization is on the order of B* y /Ep. Note that the 
factor 1/2 in B* y = I\A*\iri2k F /2 can now be interpreted 
as showing that only one half of the electron modes is 
gapped. When going to the tight-binding description this 
amplitude must in addition be weighted by the ratio of 
electron to nuclear spin densities n e i/ni, expressing that 
the number of electrons is much smaller than the num- 
ber of nuclear spins I in SWNT's and GaAs quantum 
wires. This leads to our estimate of the local electron 
polarization 



ISA*m 2 k F riei 



2E F 



^ cos(2kFri)' 
sin(2fcjr r'i) 

V 



(65) 



This argument gives us furthermore an upper bound to 
the effective hyperfine coupling constant, IA* < Ep/2, 
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telling that at most half of the electrons can be polarized. 
Further bounds and self-consistency checks are discussed 
in Sec. ED 

For the chosen systems, this electron polarization is 
small. Assuming ni2k F = 1, we find for GaAs quantum 
wires | (Sj > |/5 = 2 x 1(T 3 and for SWNTs \{8 t )\/S = 
3 x 1(T 6 (see Table U). 

E. Feedback on nuclear spins 

The partial helical polarization of the electrons natu- 
rally modifies the susceptibilities and so the RKKY in- 
teraction. In addition, according to the principle actio = 
reactio, the polarized electrons together with the renor- 
malized coupling constant A* create a magnetic field that 
acts back on the nuclear spins, hence polarizing them. 
Therefore the stabilization of the nuclear order has now 
two ingredients: The minimum of the RKKY interaction 
as before (yet with a modified shape), and the Zeeman- 
likc (but hclimagnetic) polarization by A* (Si). 

We stress that these are two different energy scales and, 
in particular, the Zeeman-like energy does not provide an 
upper bound to the RKKY scale. Indeed, both expres- 
sions are of order Aq because (S) ~ A*/Ep as shown 
just above, and so they do not follow from a first and 
a second order perturbative expansion. In addition, we 
have seen in the previous section that the electron po- 
larizations are indeed very small. Therefore we shall sec 
below that the RKKY interaction is dominated by the 
gapless (f>- modes, and hence involves different electrons 
than (S). The bounds for the validity of the perturba- 
tion theory are, in fact, imposed differently. As this is 
an important criterion of controllability of the theory, we 
analyze it in Sec. IVI Bl We show there that the per- 
turbative expansion is indeed justified for SWNTs and 
GaAs quantum wires. 

The modified Hamiltonian for the nuclear spins then 
becomes 




In the derivation of the modified RKKY interaction J' 
we suppress any occurrence of (Si) because such terms 
are of order O(Aq) and are neglected in the perturbative 
expansion. Fluctuations involving the gapped fields </>+ 
and 8 + have furthermore a much reduced amplitude due 
to the gap I A* , and can be neglected compared with the 
RKKY interaction carried by the gapless modes 0- and 
only. Fluctuations of the gapped fields become in fact 
important only at temperatures fc^T > I A*, which is 
much larger than the characteristic temperatures of the 
ordered phase (see Table IJ) . This allows us to neglect 
any occurrence of + and 6+ in the RKKY interaction. 
The details of the modification of the susceptibilities are 
worked out in Appendix IA 21 

The result is a susceptibility, and so a J g , of a gapless 
LL described by 0_ and of the same form as Eq. (|2"5|) 



with modified exponents g a — > g' a that are determined by 
the prefactors of the </>_ and 6L fields in the transforma- 
tions (j44|) -([47 |) . and a modified velocity vf — * V-. Since 
the nuclear Overhauser field singles out the spin (x,y) 
plane over the z direction, anisotropy appears between 
Xx = Xy and Xz- This is mainly expressed in different 
exponents g' x y ^ g' z , but also in that the amplitudes of 
Xx,y are only 1/2 of that of \z because one half of the 
correlators determining Xx,y depend only on + , while 
all correlators for Xz depend on 0_ and 

From the results of the detailed calculation in Ap- 
pendix |X2] we then see that the new RKKY interaction 
J' q a has precisely the same form as Eq. (f2"3"|) with the 
replacements 

J' q x > y = J*> y (g' Xty ,V-)/2, ^ = J 5 U»-), (67) 

and the exponents 

g' = g' Xty = 2K C /K S (K C + K~ l ), (68) 
g' z = {K C K~ X + K C K S )/2(K C + A'7 1 ), (69) 

satisfying g' XjVjZ < g and g' z < g' xy (for K s = 1 we have 
precisely g' z = g' xy /2) for the nanotube and quantum 
wire systems. Note that this single-band result is also 
quantitatively valid for the SWNTs, as explained in Sec. 

Em 

For the exponents we have quite generally g x y > g' z . 
Together with the difference in amplitudes we see that 
I J' q x ' v \ < \J' q z \- Naively this would mean that the system 
could gain RKKY energy by aligning the nuclear spins 
along the z axis. However, this would destroy the feed- 
back effect, and so lead to a large overall cost in energy. 
The helical order in the (x, y) plane is therefore protected 
against fluctuations in the z direction. Since this is an- 
other important self-consistency check, a detailed analy- 
sis can be found in Sec. IVI CI 



F. Modification of the cross-over temperature 

The analysis above shows that the ground state magne- 
tization of the nuclear spins remains a nuclear spin helix 
confined in the spin (x, y) plane even when the feedback 
is taken into account. The order parameter for the he- 
lical order remains rri2k F ■ Thermally excited magnons 
reduce this order parameter in the same way as before. 
As long as the gap B xy remains much larger than fc^T the 
gapped modes 4>+ and 6 + remain entirely frozen out, and 
the reaction of the electron system is described by only 
the ungapped modes 0_ and 9-, leading to the RKKY 
interaction J' q a as derived just above in Sec. IV El The 
evaluation of the magnon occupation number leading to 
Eq. (|38[) remains otherwise identical upon the replace- 
ment J 2 fe F — * J'^ F with a = x,y. Since |Jf| > |Jq' y |, 

(2) 

there is now a second magnon band uj q with negative 
energies at q f» 2kp (see Appendix [Cj. This usually 
means that the assumed ground state is unstable. But 
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in the same way as stated in the previous paragraph, the 
feedback protects the ordered state against such desta- 
bilizing fluctuations. The details are again worked out 
in Sec. IVIC1 and as a consequence we can neglect this 

(2) 

second magnon band ui q entirely. 

Let us now look at the influence of the electron po- 
larization on the nuclear spins. The effective magnetic 
field created by the electrons is B e i = \A*(Si)\. With the 
polarizations estimated in Eq. ([55)1 we obtain for GaAs 
quantum wires B e [ ~ 0.3 /.teV ~ 4 mK, and for SWNTs 
B e i ~ 36 neV ~ 0.4 fiK. Both scales are very small com- 
pared with the values of T* we shall obtain from the 
modified RKKY interaction right below. Hence, we can 
entirely neglect these magnetic fields, and so the gap they 
generate in the magnon spectrum. 

The magnon band is thus of the same type as Eq. (j34|) . 
and repeating the analysis of Sec. HVCI we obtain a mag- 
netization of the form 



3-2g' 



(70) 



with the cross-over temperature [combining Eqs. (|39|) 
and (JS7J)] 



k B T* =I\A \D' 



where g' is given by Eq. 
©], and 

l — 2g' 

\v Fj 



l-2g' 
3~2g' 



(71) 



I\A 

, again A Q = Hvp/a [Eq 



3 in(7r g ')r 2 (l~ g ')(2^) 2 g'~ 4 r 2 ( g '/2) 
2 r2(l- 5 '/2) 



(72) 

with u_ from Eq. (|5 1[) . These expressions replace the 
magnetization ra^kw m Eq. (|38[) and the temperature 
T * in Eq. (|39[) . As noted in Sec. IV Bl the m^kp can 
be directly detected by measuring the electron excitation 
gap B* y . 



VI. SELF-CONSISTENCY CONDITIONS AND 
GENERALIZATIONS 

The strong renormalization of the system properties 
through the feedback between the electron and nuclear 
systems below T* requires a reexamination of the under- 
lying conditions. We start with discussing the validity 
of the LL theory, which forms the starting point of our 
analysis. The validity of the renormalized RKKY treat- 
ment is examined next. This is followed by the inves- 
tigation of the stability of the nuclear helimagnet to a 
macroscopic realignment of the nuclear spins in the z di- 
rection that seems to be favored by the anisotropy of the 



modified RKKY. We show that the nuclear helimagnet 
is stabilized through the feedback. Finally, we show that 
intrinsic anisotropy in the hypcrfine interaction does not 
change our conclusions as long as it maintains a finite 
magnetization along a cross-section through the ID con- 
ductor. The validity of using a single-band model for 
SWNTs is discussed in Sec. IVTTll 



A. Validity of Luttinger liquid theory 

The LL theory defined by Eq. ([22|) is an exact theory 
for ID electron conductors with a perfectly linear electron 
dispersion relation. The eigenstates of such a system are 
bosonic density waves. Including electron-electron inter- 
actions does not change the nature of these eigenstates, 
but leads mainly to a renormalization of the LL param- 
eters K c and K s . Realistically, however, the dispersion 
relation is not perfectly linear and restricted to a finite 
bandwidth, and electron-electron interactions can have a 
more substantial influence. In such a situation the LL 
theory remains valid as long as (for the considered physi- 
cal quantities) the bosonic density waves remain close to 
the true eigenstates and decay only over a length scale 
exceeding the system length. 

Deviations from LL behavior induced by the electron 
band curvature close to ±kp was investigated in Refs. 
I86ll87ll88l . Let us encode this curvature in a mass m* such 
that the electron dispersion reads e q = Vpq + q 2 /2m* , 
where q is measured from ±fcp- Defining then the pa- 
rameter e = {w — e q )/(q 2 /2m*) it was shown^ that de- 
viations from LL behavior become important at |e| < 1 
provided that |g| <C kp. In our case, the electron correla- 
tion functions are evaluated in the static limit tu = 0, and 
by setting q m = 2vpm* this condition becomes |q| > q m . 

For armchair SWNTs we estimate^ 9 - m* « 0.2rao (with 
mo the bare electron mass) within a few 0.1 cV about 
the Dirac points. This mass is very large, reflecting the 
almost perfect linear dispersion of the armchair SWNTs. 
Accordingly this leads to a q m « 3 nm _1 3> hp, or to an 
energy scale of about 2 eV. Hence, q m ^> \q\ for any \q\ <C 
kp. For GaAs quantum wires the effective mass at the 
T point is2£ m* = 0.067m , and so q m = 2 x m vp/h = 
2.3 x 10 8 m _1 . This is slightly larger than kp and again 
we find that q m ^> \q\ for any \q\ -c kp. Therefore, for 
both systems the curvature-induced deviations from the 
LL theory are negligible. 

A different curvature-induced deviation from standard 
LL theory occurs at very low electron densities, leading 
to the so-called incoherent LL (see e.g. Refs. l58ll9lU92ll93l 
and references therein). At these densities the Coulomb 
energy E pot largely overrules the kinetic energy E km of 
the electrons (expressed by a ratio R = E pot /E kln 3> 1), 
and the electrons order in a Wigner crystal, while the 
electron spins form a Heisenberg chain. Such a sys- 
tem still allows a bosonized description , 91 ' 92 yet with a 
large splitting between spin and charge excitation en- 
ergies. Realistic temperatures lie above the spin ex- 
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citation energies, but can lie below the charge excita- 
tion energies. The charge fluctuations then remains in 
a LL state, while the spin fluctuations have an incoher- 
ent behavior. The ratio R depends much on the band 
mass m of the system. For a quadratic dispersion we 
have E kln = k 2 F /2m oc n 2 el /2m. For a potential energy 
E pot = n e ie 2 /e (with e the electron charge and e the 
dielectric constant) we obtain R — E pot / E km oc m/n e i. 
The incoherent LL regime can therefore be reached in sys- 
tems with a large mass and a low electron density^ In 
the GaAs/AlGaAs heterostructure of the quantum wires 
we have2£ e « 12, and E kln = Ep. With the values from 
Table |T] we find that R = n e ie 2 /eEp ^0.2, excluding the 
incoherent LL. 

For a linear spectrum such as in the SWNTs the cri- 
terion above does not apply. Indeed, for a linear dis- 
persion E kln oc n e i and so the ratio R = E pot /E kln is 
independent of density. Electron interactions then pri- 
marily modify the dielectric constant e and so K c . This 
leads to a much weaker dependence of R on the inter- 
action strength, and makes the LL description valid for, 
in principle, arbitrary electron-electron interactions. For 
SWNTs we have the estimate^^ e 1.4 leading to 
R, = e 2 /irhvpe « 0.6, allowing us to exclude the inco- 
herent LL as well. 



B. Validity of RKKY approximation and bounds 

on T* 

The RKKY approximation is valid under two condi- 
tions. First, as it is a perturbativc expansion in powers 
of Aq, we must verify that higher perturbative orders 
remain smaller than the lower perturbative orders. Re- 
lated to this we must examine that the energy scale k B T* 
obtained from the RKKY interaction does not violate 
bounds imposed by the original Hamiltonian ([7]). Sec- 
ond, the separation of time scales between the electron 
and nuclear spin systems must be guaranteed, in order 
to be allowed to interpret the RKKY interaction J' q as 
instantaneous for the nuclear spins. 



1. Upper bound on T* 

The perturbative derivation of the RKKY 
interaction^ or equivalently its derivation through 
a Schricffcr- Wolff transformation^ consists in an ex- 
pansion in Aq. The lowest order is proportional to Aq, 
the scale of the second order is set by J q oc A^/Ep, 
higher orders scale in further powers of Ag/Ep, and 
so the condition of validity of perturbation theory 
is usually set equal to the condition |Ao|/-Ef <C 1. 
This condition is perfectly met for the SWNTs or 
the GaAs quantum wires. Yet, there are a few sub- 
tleties. First, in the absence of the feedback, we have 
J q = (Aq/Ep)x (number), and the latter number can 
become very large. Its maximum defines indeed the scale 



UbTq. Through the Schricffcr- Wolff transformation we 
have in addition eliminated the term in the Hamiltonian, 
which is linear in Aq, and so there is no longer a proper 
"first order" expression to which we can compare J q . 
Hence, the validity of the RKKY scale fc^Tg must be 
checked in a different way. 

If we entirely neglect the electron Hamiltonian, i.e. 
H = AqS{ - li, we can ask which maximal energy scale 
can be obtained from the hyperfine Hamiltonian for the 
nuclear spins. Obviously this scale is obtained by polar- 
izing all electrons such that H = J2iAo{n e i/ni)Se ■ Ij, 
with 5* = 1/2, n e i/ni the ratio of electron to nuclear 
spin densities, and e an arbitrary unit vector (that may 
or may not be position-dependent). For temperatures 
smaller than the field = \ Ao\Sn e i/ni all nuclear spins 
are aligned along e because a mismatch with the fully 
polarized electrons costs the on-site energy B'q. Essen- 
tial for this argument is that flipping a nuclear spin out 
of its alignment costs only energy from the hyperfine in- 
teraction. The electron system is assumed to be energeti- 
cally unaffected by this process (the hyperfine interaction 
conserves the total of nuclear and electron spins, and so 
the electron spin changes as well). Otherwise said, the 
electron Hamiltonian H e i is independent of the electron 
polarization. Within this framework the RKKY coupling 
between the nuclear and electron spins corresponds just 
to a more sophisticated way of treating the nuclear spin 
fluctuations. Since the electron state has no influence by 
assumption, the maximal energy scale set by Bq cannot 
be overcome by any characteristic temperature obtained 
through the RKKY interaction. This argument therefore 
applies directly to the case when we neglect the feedback 
between the nuclear spins and the LL. The condition then 
becomes 

k B T* <B*=SA ^. (73) 
m 

With the values from Table [T] wc can verify that this 
condition is indeed satisfied. 

The situation is of course very different if both sys- 
tems are tightly bound together. In this case, one pays 
not only the energy Bq but also the energy resulting from 
the modification of the electron state. Through the feed- 
back between both systems, this extra cost in energy is 
roughly taken into account through the renormalized hy- 
perfine coupling constant A* [see Eq. (|57|) ]. If wc assume 
that A* fully describes the maximal electron response to 
the hyperfine coupling, then the scales obtained from the 
modified RKKY description again cannot overcome this 
scale. Hence, we have the modified condition 

77 / 

k B T* < B* = SA* — . (74) 
ni 

We must interpret this inequality with caution though. 
Only one half of the low-energy electron modes contribute 
to the renormalization of A* . The Coulomb interaction 
between the remaining electrons, which in fact substan- 
tially modifies the RKKY interaction and determines T* , 
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is not taken into account. Yet, this modified RKKY in- 
teraction is a direct consequence of the strong coupling 
to the nuclear system as well, and so the scale B* may 
still require further adjustment. A hint for this is seen for 
instance in Fig. [6l where k B T* exceeds B* for K c < 0.5. 

Nonetheless we use Eq. (|74]) as an upper bound to 
ksT* , since we do not know if such an extrapolation 
beyond B* remains valid within the RKKY framework. 
However, we interpret Eq. (|74|) as assuring the validity 
of the theory, and not necessarily as a maximal upper 
bound on k B T* . 

If we look again at Table |T] we see that k B T* w B* for 
the selected values of GaAs quantum wires and SWNTs 
[within the 0(1) uncertainty]. Hence, the fluctuations 
of the gapless modes stabilize the nuclear order up to 
the scale B* set by the gapped modes. This equality 
of the scales for the parameters of Table Q] is actually a 
coincidence, as can be seen from Fig. |6] 

Note that if B* is controlled by the cutoff £ = L, 
increasing the system length also increases B* because 
more electrons are involved in the feedback effect. In 
SWNTs, for instance, £ = L and B* cx (i/a) 1 ^ = 
(L/a) 0A . This means that doubling L corresponds to 
an increase of B* by about 1.3. Since T* is independent 
of L such a control of the bound B* may be quite useful. 

2. Separation of time scales 

The RKKY interaction requires that the electron re- 
sponse to a change of the nuclear spin configuration is 
instantaneous, and so a strict separation of time scales 
between both systems is mandatory. The dynamics of 
the nuclear spins is described by the Hamiltonian ([66)) 
and consists of two parts. The precession of the nuclear 
spins in the magnetic field generated (sclf-consistently) 
by the polarized electrons and the renormalized cou- 
pling constant A*, and the dynamics from the RKKY 
interaction J'. The former leads to an energy scale 
A*(S) ~ (A*) 2 /E F , which needs to be compared with 
Ep. This results in the condition 

A* < E F . (75) 

Notice that this argument is very different from the pre- 
vious argument leading to Eqs. (|73|) and ([74]) . as it re- 
quires the physical, fully self-consistent averages, not a 
maximal energy condition. From Table[T]we see that con- 
dition (fTS")) is met for SWNTs and GaAs quantum wires. 
Since A* /Ep measures essentially the proportion of elec- 
tron spin polarization (see Sec. IV Dp . it means that only 
a very small fraction of the electrons is polarized. 

On the other hand, the time scale set by J' q can be 
identified with the dynamics of the fluctuations it de- 
scribes, and so with the magnon dynamics. We therefore 
compare the maximal magnon velocity with the Fermi 
velocity vf- For temperatures T < T*, the maximal 
magnon velocity v m is obtained by the slope of u q at the 
momentum cutoff q = tt/L. Wc can then use the T = 



expression for J q , as it sets an upper bound to the slope, 
and obtain from Eqs. ([30]). ([34]) . and j67|) 

vjn ^ (2-2 g ')sin(^')(2/7r) 3 - 2 gV I A% f L\ 3 ~ 29 ' 

(76) 

The first 3 factors are small and can overcome the large 
last factor. Indeed, for SWNTs we have v m jvp ~ 
10 _9 /iVj_, and for GaAs quantum wires (with L = 10 
/xm) Vm/vp ~ 10~ 3 /N±. The already small prefactors 
are in addition strongly suppressed by the number N± of 
nuclear spins in the direction across the ID system. The 
small values for v m mean that q = tt/L lies already in a 
region where uj q is essentially flat. This is, in fact, the 
same criterion we have used for the determination of T*. 

The necessary separation of time scales is therefore ful- 
filled by the systems under consideration. 

C. Stability of the planar magnetic order 

We have observed in Sec. IV El that quite generally 
\J' q z \ < \J' q x ' v \ because of the smaller exponent g' z and 
the larger prefactor of J' q z . Naively, this means that the 
nuclear spin system can gain RKKY energy by forming 
an Ising-like configuration along the spin z direction. An 
alignment in this direction, however, would destroy the 
feedback described above, and so destroy the net energy 
gain from the planar (x, y) order in both the electron and 
spin systems. To keep this feedback and the energy gain 
active, a deviation from the planar order of the nuclear 
spins is not possible. 

Indeed, the narrow minimum of J' q at q — 2kp im- 
plies that if there is a magnetization m z along the spin 
z direction it has only q = ±2fci? Fourier components, 
because they are energetically most favorable. The real- 
ity of the expectation value of each nuclear spin then 
imposes that we can write the ground state expecta- 
tion values in the form = IN±m2k F [cos(2kpr)e x + 
s'm(2kpr)e y ]+IN±m z sm(2kpr+rj)e z , where < m z < 1 
is the magnetization along the z direction and rj an arbi- 
trary phase. For m z ^ we see that a full polarization 
= IN± is no longer possible in general. Instead, the 
maximally possible polarization is determined by the con- 
dition m\ kF + m 2 , = 1. Choosing such a rri2k F also mini- 
mizes the nuclear spin energy for a fixed m z , so that we 
can consider the latter condition as being fulfilled when 
seeking the absolute ground state energy. 

Using the bosonization approach^ in the same way as 
in Sec. IV El we see that the new m z component leads 
to an Ovcrhauscr field for the electron system of the 
form B z sin(v / 20 c ) sin(v / 20s) with B z = IAom z which 
is relevant under the RG, plus oscillating terms which 
can be neglected. This new term competes against the 
term B xy cos(\/2K (f> + ) = B xy cos(v / 2((/) c + S )) (with 
B xy = IA iri2k F ) hi the RG, because it involves the field 
(j) s which is conjugate to 8 S . By the uncertainty principle, 
the pinning of both fields is impossible, and generally the 
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term with the larger amplitude, B xy or B z , dominates the 
RG flow to the strong coupling fixed point 

To our knowledge there is no method allowing a pre- 
cise evaluation of this RG flow. The following esti- 
mate, however, is sufficient for our needs. At the cut- 
off scale determining B* y , both B xy and B z have flown 
to strong coupling, although through the competition 
we have B* y (m z ^ 0) < B* y (m z = 0). The scal- 
ing dimension of B z is 2 — K c /2 — K s /2 and so for 
K s = 1 identical to the scaling dimension of B xy [Eq. 
(|52|) ]. Hence, the initial ratio B z /B xy — m z /m2k F re- 
mains, up to small corrections, constant throughout the 
RG flow. This allows us to estimate the decrease of 
B* y to be proportional to m z , which leads to a cost 
in the electron ground state energy per lattice site on 
the order of A£ cost = B* y - B* y (m z ) ~ m z B* xy [with 
B* y ~ B* y (m z = 0)]. This cost must be compared 
with the gain in nuclear spin energy per lattice site, 
given by AE^ ~ I 2 m^l F - & ~ ^"MU 
Since A#s ain - m 2 z and AE cost ~ m z the relation 
A£s ain < A£ cost is most likely for small m z and defi- 
nitely always fulfilled if I 2 \J2k F \ < B xy . Since close to 
T* [see Table H] we have B* y > k B T* ~ I 2 \J' 2 l F l we 
conclude that any m z ^ is energetically highly unfa- 
vorable and so the helical order in the spin (x, y) plane 
is stable. 

With this argument we also see that the second 
magnon band derived in Eq. (fC7j> . ) = 2I(J z q - 

(2) 

J%k F )/N±, which has io q < at q ~ 2fcp, is of no im- 
portance. Any macroscopic occupation of these negative 
energy states, which would normally signify an instability 
of the assumed ground state, is energetically forbidden. 

(2) 

The remaining u> q > describe 2kp fluctuations in the 
z direction. Their effect is the same as the fluctuations in 
the (x, y) plane described by the first magnon band ujq 1 ^ 
in Eq. (|C6[) . yet they involve large momenta ~ 2kp. Due 
to this we can neglect this second magnon band entirely 
in this theory. 

Let us finally note that a pure m z magnetization (with 
the (x, y) component rri2h F = 0) would lead to a similar 
feedback effect as the helical magnetization and open up 
a gap in the <p c + 4> s sector. This would lead to a spa- 
tially oscillating Ising-like average magnetization in the 
z direction. In contrast to the (x, y) helical magnetiza- 
tion, the condition = NI is then fulfilled only when 
cos{2kpTi) = 1. The resulting nuclear magnetic energy 
lies therefore much above the energy from the (x, y) mag- 
netization; the minimum is only about half as deep as in 
the latter case, and thus such a state is not assumed by 
the system in the ground state. 



D. Anisotropic hyperfine coupling 

The strong feedback between the nuclear spins and the 
electrons occurs only if there is a nonzero Ovcrhauscr 
field on every cross-section through the ID conductor. 




FIG. 10: Illustration of the nuclear spin configuration min- 
imizing the energy of Eq. (|80|) for the case of a complete 
electron polarization S pointing upwards in the figure. The 
figure shows how the sign of the hyperfine constant changes 
when going around the nanotube cross-section. The overall 
magnetization of this configuration is m = —0.17 (along x, 
normalized to —1 < m < 1). See also Ref. l28l . 



This was ensured by the ferromagnetic locking of the nu- 
clear spins by the coupling to the single transverse elec- 
tron mode. Anisotropy in the hyperfine interaction can 
perturb this situation, and a reinvestigation of this very 
important first assumption on the nuclear spins becomes 
mandatory. 

Here we focus on the case of carbon nanotubes, in 
which indeed anisotropy is present through the dipolar 
interaction between the electron and nuclear spins on 
the curved surfaced Rotational symmetry imposes that 
if anisotropy is present, it occurs between the radial di- 
rection (r) , the tangential direction on the circular cross 
section (t), and the direction along the tube axis (c). 
Writing the nuclear and electron spin operators in this 
local basis, I, = and S, = (Sl,Sj,Sf), the 

hyperfine Hamiltonian can be written as 



H = J2 [A r IlSl + AHtSt + A% c Sf] , (77) 



where i runs over the 3D nuclear spin lattice. We as- 
sume henceforth thatSS A r = -2A* = -2A C = 2A Q > 0, 
and neglect the A c term as it turns out to be smaller than 
the couplings in the plane spanned by the r and t compo- 
nents. Let us identify this plane with the spin (x, y) plane 
and rewrite the local component in an global spin basis as 
ST = Sf cos(C i )+af 8in(Ci), S* t = -Sf sin(G)+aj'cos(C i ) ! 
and analogously for the Ii operators, where Q is the polar 
angle of site i on the circular cross-section (see Fig. [T0| . 
The Hamiltonian can then be written in the form 
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H = A Q 



^(Sf, Sf) \l + 3cos(2Ci)cr z +3siii(2C0<r a 



(78) 



with er Xj2 the Pauli matrices. As in Sec. IIII Bl we as- 
sume that the electrons are confined in a single transverse 
mode, which allows us to sum over the transverse compo- 
nents of the indices i for the nuclear spins. We see then 
that the electron spin couples only to the ferromagnetic 
component of the total nuclear spin on the circle, Iff. = 

= E i± 008(20)/^ 



E^ )V and the 2Q modes I, 
and/",, = EiX2^, V 



Since S: 



we have 
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ja 



(79) 

The coupling to the electron spin, therefore, acts simulta- 
neously on these three nuclear spin modes only. For the 
feedback it is essential that these modes carry a finite 
magnetization. This is actually the case as we can see by 
assuming an electron spin at a cross section i|| polarized 
in the S x = +1/2 direction. The Hamiltonian for 
the nuclear spins on this cross-section then becomes 



AoS 

N 



— V(V T y H 1 +3eos(2£i x ) ] 
I 2-< \ 'U'W ^ 3sin(2C 4 J J 



(80) 

The energy is minimized when at each position (i\\,ix) 
the nuclear spin is opposite to the vector on the far right 
in Eq. ([8H)l . This configuration is shown in Fig. |TD] 
and has a net magnetization along the x direction of 
m = —0.17. Hence, the coupling to the single trans- 
verse electron mode enforces a nuclear Overhauser field 
of the same type as in the isotropic case. 

Through the reduced magnetic field, we see from Eq. 
([55]) that the gap B* y is smaller compared with the 
isotropic value by (0.17) 1/(2 " s) w 0.2. This affects the 
feedback only if the smaller gap becomes comparable to 
k B T*. Since, however, B* ~ 10 fxeV ~ 100 mK for the 
isotropic case (see Table Qj, even with this reduction, it 
remains on a scale that is larger than T* ~ 10 mK. Yet it 
also lowers the limiting scale B* below k B T*, and so the 
true T* may be a bit lower as well. However, since B* de- 
pends directly on L for SWNTs [see Sec. I VI B| . choosing 
a longer sample length will push B* up to higher values, 
and so T* may keep its original value. 



VII. RENORMALIZATION ABOVE T* 

For temperatures above T* thermal fluctuations de- 
stroy the nuclear helimagnet, but there remains the sepa- 
ration of time scales between the nuclear and electron sys- 
tems. This implies that any random fluctuation at time t 



into a nonzero Fourier component rri2k F (t) of the nuclear 
magnetization triggers a renormalization of the electronic 
properties in the same way as before and leads to an in- 
stantaneous Overhauser field B* y (t) = Im 2 k F {t)A* (t)/2. 

Any detection of this Overhauser field would have to 
rely on the measurement of nonvanishing quantities such 
as (B* y ) 2 . We emphasize though that the renormaliza- 
tion occurs only for the q = ±2fcf Fourier modes of the 
magnetization. Close to T* thermal fluctuations may 
still occur preferably in the vicinity of the minimum of 
J q at q — ±2fcp. However, as T is raised away from T* 
the relative weight of the q = ±2kp modes with respect 
to the other (N — 2) Fourier modes rapidly drops and 
approaches \/N at high temperatures. Except for a T 
in the close vicinity of T* the renormalized Overhauser 
field is probably not detectable. 

On the other hand, as long as k B T remains below the 
maximally possible B* y = I\A*\/2, the m 2 k F fluctua- 
tions sporadically open a gap for one half of the electron 
modes whenever B* y (t) exceeds k B T. This temporar- 
ily freezes out of one half of the conduction channels, 7 - 
and so reduces the time averaged electric conductance 
of the system. The reduction is given by a temperature 
dependent factor / between 1/2 and 1, with the limits 
/ = 1/2 at T — > T* and / = 1 when crossing through 
k B T - I\A*\/2. Note that if the correlation length £ is 
given by the system length L, increasing L also increases 
A* [see Eq. ([57])]. and so the onset of / < 1 occurs al- 
ready at higher temperatures. 



VIII. VALIDITY OF THE ONE-BAND 
DESCRIPTION FOR CARBON NANOTUBES 

The band structure of armchair carbon nanotubes con- 
sists of two Dirac cones (or two valleys) centered at the 
momenta^! K = ±2ir/3a (folded into the first Brillouin 
zone). Hence, the LL theory depends not only on spin 
and L or R movers but also on an index a — 1, 2 labeling 
the two Dirac cones. We have neglected this additional 
quantum number a in the previous sections with the ar- 
gument that, once the feedback is included, the results 
obtained from the single band (single cone) description 
employed so far are quantitatively the same. In this sec- 
tion this shall be explicitly shown. 



A. Luttinger liquid theory for nanotubes 

Without electron-electron interactions SWNTs of 
the armchair type have the effective low energy 
HamiltonianSS 

(81) 

where ipiaa is the electron operator of I = L, R movers 
in cone a with spin tr, and i/j a — tytaa ' 1S the full 
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electron operator. This theory can be bosonized within 
each cone in the usual wa y 68 i 69 ' 70 by introducing the 
boson fields <f> va and ua with v = c, s. In particular 
-V0 CQ \/2/7r = p( a , where p{ a = T, ia ^L^^a is the 
forward-scattering part of the density operator in cone 
a. The basic bosonization identity (|20p is now enriched 
by the index a: 



4>laa 



tya 



i[ek F +ak v ]r+-if[<p c . a -e6 ca +a( 



-te. a )] 



(82) 

where kp is measured from the crossing point of the two 
branches of the corresponding Dirac cone, and so the 
position of these Dirac cones in the Brillouin zone has 
to be taken into account through the additional phase 
factors ±ik v r with k v = 47r/3a (modulo any reciprocal 
lattice vector). 

The electron-electron interactions considered here are 
density-density (Coulomb) interactions and can be split 
into a part depending only on the forward-scattering 
densities p{ a and a term depending on backscattering 
I — > —£, inter-cone scattering a — ► —a, and Umk- 
lapp scattering. The latter processes lead to a gap in 
the system, but with a magnitude that is exponentially 
suppressed with an increasing diameter of the tube^ 
At realistic temperatures this gap is only important for 
very narrow tubes and can be excluded here. Indeed, 
from Ref. [7^ the dominating gap is determined by 
777.5 = wq exp(— nhvp / ^/2b), where too = 7.4 eV is the 
bandwidth of the it electrons, and b sa 0.1 ae 2 / R is 
the inter-valley scattering amplitude with e the electron 
charge and R the radius of the SWNT. For (n, n) arm- 
chair SWNTs we have^ R = | C|/2tt = V3an/2n, with 
C the chiral vector (see also Sec. IIIIBI) . Putting in num- 
bers (see Table |I| we find rrif, ~ u>q exp(— 2.2n), and for 
n = 13 we have mb ~ 3 peV ~ 30 nK. We notice, how- 
ever, that b has an order 1 uncertainty^ and can become 
larger for very well screened interactions. LL physics at 
the millikelvin temperatures considered here is valid for 
gaps that lie below these temperatures, which therefore 
requires a sufficiently wide nanotube and a not too short 
screening length (yet short enough to allow a local de- 
scription for the interaction). In addition, the experi- 
ments of Ref. H3 show no indication for such a gap down 
to a temperature T = 1.6 K, confirming the LL picture. 

The system is then describable by a LL where only 
the forward scattering interaction part remains of im- 
portance. It is of the form 

V = J drdr'V{r - r')[p f cl {r) + p f c2 (r)}[p f cl (r') + p{ 2 {r')\ 



drC/[V0 cl (r)+V^ c2 (r)] 5 



(83) 



where U can be related to a charging energy and V(r) is 
sufficiently screened to allow a local description (yet see 
the preceding paragraph). 

The interacting Hamiltonian can then be diagonal- 
ized through boson fields that describe the symmetric 



S and antisymmetric A parts of density fluctuations be- 
tween the a = 1,2 cones, <f) v s = {4>ui + 4>v2)/\/% and 
4>vA = (0i/i — 0i^2)/v / 2- Similar definitions hold for the 
9 ua fields. The electron-electron interaction acts only on 
the symmetric (c, S) sector. Therefore the LL theory re- 
mains that of a noninteracting system in the (c, A), (s, S), 
and (s,A) sectors (more precisely, it remains that of al- 
most a noninteracting system due to slight renormal- 
izations through the backscattering interaction). The 
bosonic Hamiltonian for the nanotube then takes the 
form 



(v^p) 2 + v v pK vP (ve vP f 



(84) 

where P = S, A. Here K c a ~ K s g m K s a ~ 1, while K c s 
is strongly renormalizcd by the interaction and typically 
is abou t 68 ' 69 i 70 K c s ~ 0.2. The velocities are given by 
v vP = v F /K vP . 

Correlation functions computed with this theory are 
consequently described by power laws with exponents 
composed of all four LL parameters K u p. Of particu- 
lar interest for us are the backscattering amplitudes ap- 
pearing in the spin operators S x and S v , for instance for 
spin-flip backscattering R — > L, 



_ e — ■ e 
for inter-cone scattering a = 1 — ► 2, 

4 2 T^iU ~ e 2lfc " r e75 (0t2T - 0ili) 

_ e 2ik v r e i(-<t>cA-0cA+$sS+9ss) 
V'LtV'RU ~ e 2ifc„r e ^|(teT-te|) 

_ e 2ik v r e -i(-<t>cA+6eA+<f>,s-0,s) 



^2ik F r a -ys(4>La,f+<pRc,i) 



2ik F r i{<t> s+a4>oA+SsS+a.B sA ) 



e 2ik F r e iV2(4, ca +9 3a ) 



(85) 



(86) 



(87) 



and for the combinations L — > R together with a = 1 
2 or a = 2 -> 1, 



V'Lj.V'iiiT 



3 2i(fc F +fc„)r 



e V2' 



e 2 l (fc J? + fc„)r e j( 



*i2T+0-Rli) 
;-e cA -4, sA +8 sS ) 



e 2i(k F +k v )r e ^(<l>L2i+<t>Rir) 

e 2i(k F +k v )r e i(<p a s-8 cA +<p sA -e sS ) 

e 2i(k F -k l ,)r c ^(4>LlT+<pR2i) 



e 2z(k F -k v )r e i(<p c s+8 cA +<p sA +e sS ) 



^2i{k F —k v )r 



e 2i(k F -k v )r e i(<l>cS+6cA-<t>zA-8ss) 



(89) 

(90) 
(91) 



All these expressions enter the susceptibility and lead 
to power-law divergences at the momenta 2kp, 2k v and 
2(fcp ik v ). The operators involving backscattering L 
R depend all on <f> c s and 3 of the other fields, leading to 
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the same exponents g — (K c s + 3)/4 0.8. Eqs. (|86|) 
and (|57|) do not depend on 4> c s and so lead to correlators 
with a larger exponent, and so a shallower minimum for 
the RKKY interaction. 

This means that the RKKY interaction J q between 
the nuclear spins has 3 equal minima and so the helical 
order is in principle not a simple hclimagnct, but involves 
the superposition of the 3 spatial frequencies. Yet this, 
once again, neglects the feedback on the electron system, 
which selects only the 2kp frequency as relevant for the 
feedback and for lowering the ground state energy, as we 
shall see next. 



B. Feedback from the nuclear Overhauser field 

Let us assume that the Overhauser field consists of the 
spatial frequencies 2kp and 2(kp ± k v ) associated with 
the minima of the susceptibility. For completeness we 
keep also here the local minimum at 2k v , and so write 
the coupling of the nuclear magnetic field to the electron 
spin as 



different for the B$ and B4 processes. Because of the 
mixing between the Dirac cones, the processes L ^<-> R j 
and L l*-* R 1 have the different momentum transfers 
of 2{kp + k v ) and 2(kp — k v ), respectively. The open- 
ing of a gap for one of those processes entirely eliminates 
the steep minimum of Xx.y{q) at the corresponding mo- 
mentum q = 2(kp ± k v ), and so eliminates this energy 
minimum of the nuclear spins. Consequently, the nuclear 
spins will reorient themselves into one of the remaining 
energy minima, i.e. the amplitudes B3 and B4 vanish. 
The only remaining minima are then those of B±. With 
vanishing B3 and B4, however, the corresponding gaps 
vanish as well. Hence, even though the B34 energy min- 
ima exist, they cannot be populated by the nuclear spin 
modes l q because this immediately opens gaps and costs 
energy. This is very similar to the stabilization effect 
discussed in Sec. IVI CI 

Therefore, any nuclear spin order with 2(kp ± k v ) is 
unstable, and the only stable minimum for the feedback 
is the nuclear helical order at 2kp. 



i=l 



Bi [S x {r) cos(2hr + 77,) + S y {r) sin(2fc l r + 



(92) 

with fci = kp, k 2 = k v , k 3 = kp + k v , and k^ = kp — k v , 
the corresponding amplitudes Bi oc mfc 4 proportional to 
the ki components of the magnetization, and r\\ , . . . , 774 
arbitrary angles, which we can set to zero because they 
can be absorbed by suitable shifts of the boson fields. 
Using Eqs. (|55|) - (|89"|) to express S x ' y in terms of the 
boson fields, we are led to the interaction 



B x cos(\/2(0 c i + 6 sl )) + cos(\/2(</>c2 ■ 
+ 2B 2 cos((/) cA - <j) sS ) sm(6 cA - sS ) 

+ 2B 3 COS(> c s - 9 cA ) COs((/) sA - Oss) 

+ 2B 4 cos(<p cS + c a) cos{(f> sA + 9 sS ) 
+ oscillating terms. 



"s2 



)) 



(93) 



The terms proportional to B\ are precisely the combina- 
tions + oc 0c + 9 S used previously, and lead to a gap in 
each Dirac cone separately. The terms proportional to 
B2 involve fields that are conjugate to each other, and 
so this interaction remains always critical and can be ne- 
glected. The arguments proportional to B3 and B4 are 
relevant and open gaps in the corresponding boson fields. 
Their effect is, however, quite different from that of the 
B\ part. 

Indeed, for the B\ part, the opening of the gap in 
the <p + oc cj) c + S channels leaves still the combination 
(\> c — 9 S gapless [see Eq. (|42p]. Both combinations ap- 
pear in the susceptibilities \x,y and correspond to the 
processes L R { and L R f. Hence, even with 
the gap for <fi c + 8 S there remains (at T = 0) a power-law 
singularity at 2kp in Xx,y due to (f> c — 8 S with the conse- 
quences discussed in the previous sections. This is much 



C. Validity of the single-band model 



Since the cos(v / 2(0 C a + sa )) interaction terms in each 
Dirac cone a — 1, 2 are highly RG- relevant, we can eval- 
uate their effect in each cone separately. Transforming 
into the a = 1, 2 basis, the Hamiltonian takes the form, 



I. ' a > 



+ — cos [V2((j) ca + 9 sa )j + v F {K-£ - l)V<£ c iV<A 



c2 , 

(94) 



for 



c, g, a = 1,2, and with K ca = 
K cS = 0.2, K sa = 1, and v va = 



vp/K va . Without the last Dirac cone mixing term this 
Hamiltonian is identical to two copies of the single-band 
Hamiltonian (|49|) . The coupling between the Dirac cones 
is, however, marginal in the same sense as the gradient 
products in the last term of Eq. (|49|) and has only a tiny 
influence on the opening of the gap by the cosine term 
in each cone separately. Hence, for the evaluation of the 
gap we can safely neglect this coupling. The size of the 
gap is then determined by replacing K c and K s in Sec. 
ELS by K c = K cl = K c2 w 0.2 and K s = 1. 

The resulting low-energy Hamiltonian then depends 
only on the gapless fields <f>-^ a and 9— a , defined as in 
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Eqs. (gll) 
H_ 



47j) . and is given by 



+ WF (^ 2 -1)-^V0_ 4 V^_, 2 



27V^-, 1 V0_,2 



(95) 



with 7 = «f[A7 s 2 - l]A c /2v-A s A = (1 - A 2 S )/2(1 + 
2A' 2 S ) sa 0.44. This Hamiltonian is diagonalized by 
the symmetric and antisymmetric combinations 4>s — 
(<f)- t i + (f>- :2 )/V2 and <\>a = {</>-. l - <t>-,i)/^, and simi- 
larly for the 9 fields. We obtain 



dr x - 
27 ^ 



P=S,A 



|^(V0p) 2 
K P 



v P Kp(V9p) 



1.20, A A = VT^7 



(96) 
0.75 and 



with K$ = vl + 7 
up = v-/Kp. 

The single-band model used in the previous sections 
can be applied if it produces the same results as the 
model (|96p . The essential quantities are the susceptibili- 
ties Xx,y,z evaluated with the remaining gapless fields for 
the 2kp backscattcring of Eq. (|85|) . In the single-band 
model the LL constant is A_ = 1 . From Eq. we see 
that only the combinations <ps ± 4>a and 9s ± 9a appear 
(the prefactors are the same as in the single-band model 
up to negligible corrections). The (f>s ± <\>a lead to expo- 
nents depending on (Ks + Ka)/2 f=a 0.97 and the 9s ± 9a 
lead to exponents depending on (Kg 1 + K^)/2 « 1.08. 
Since both values are very close to 1, the 2-band theory 
leads indeed to the same conclusions as the single-band 
theory. 

In fact, the single-band theory can be interpreted as 
neglecting the coupling between the Dirac cones, which is 
the approximation used in Ref. [^. As we have seen above, 
with the feedback this becomes very accurate because 
the relative coupling strength 7 between the two cones 
[Eq. (f9"B")) ] is much reduced. Indeed, corrections to the 
decoupled systems appear only at 7 2 (from expanding 
K s + Ka or Kg 1 + A" 1 ). The value of 7 w 0.44 is small 
enough compared with the original (A^ 2 — 1) w 24 in 
Eq. ([94]) so that the coupling between the Dirac cones 
has a negligible effect. 



IX. CONCLUSIONS 



the electrons and nuclear spins stabilizes this combined 
order up to temperatures that are within experimental 
reach, even though the hypcrfmc interaction generally is 
very weak. The feedback is a direct consequence of Lut- 
tinger liquid behavior and is absent for noninteracting 
electrons or Fermi liquids. 

This leads to several remarkable effects that should be 
detectable experimentally and that may be used for a 
direct proof of Luttinger liquid physics: (i) The helical 
magnetization rri2k F resulting from the nuclear spin or- 
dering follows the modified Bloch law of Eq. @. (ii) 
The helical electron spin density wave resulting from the 
renormalization triggered by the nuclear Overhauser field 
has an excitation gap proportional to rri2k F , or to some 
power of m-ikf depending on the Luttinger liquid param- 
eters [see Eq. (|59p ]. Measuring this electronic excitation 
gap is a direct way of measuring m2k F - (hi) The pin- 
ning of one half of the electron conduction modes causes 
a reduction of the electric conductance by the factor 2 
when cooling down through the cross-over temperature 
T* . (iv) Finally, the strong binding of the nuclear heli- 
magnet to the electron modes leads to anisotropy in the 
electron spin susceptibility between the (x, y) plane of 
the nuclear helimagnet and the orthogonal z direction. 

We refer the reader to Sec. [Til for a complete sum- 
mary of our results and the main conclusions. The phys- 
ical mechanism for magnetic ordering described here was 
worked out for two experimentally available systems, 13 C 
single wall nanotubes, and GaAs-based quantum wires. 
We note, however, that we expect similar physics for a 
large class of Kondo-lattice-type systems coupled to a 
Luttinger liquid defined by the conditions listed in Sec. 

IE 



Acknowledgments 

We thank B. I. Halperin and M. J. Schmidt for helpful 
discussions and critical remarks, and A. Yacoby and G. 
Barak for providing the values characterizing the GaAs 
quantum wires. This work was supported by the Swiss 
NSF, NCCR Nanoscience (Basel), and DARPA QUEST. 



APPENDIX A: ELECTRON SPIN 
SUSCEPTIBILITY 

In this appendix we evaluate the electron spin suscepti- 
bilities for the cases without and with the feedback from 
the nuclear Overhauser field. 



We have shown in this paper that the hyperfine inter- 
action between a lattice of nuclear spins and Luttinger 
liquid leads to order in both systems in the form of a 
nuclear helimagnet and a helical spin density wave for 
half of the electron modes. A strong feedback between 



1. Without feedback from Overhauser field 

Without the Overhauser field the electron system is 
described by a standard LL. 
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Wc define the static electron spin susceptibility in the 
ID tight-binding basis as [see Eq. (|T5|) ] 



a/3 
Xij 



^e-^([Sf(i),Sf(0)]), (Al) 



with an infinitesimal r\ > 0. For the further calculation it 
is more convenient to pass to the continuum description 
where the tight binding operators Sf{t)/a are replaced 
by the fields S a (r,t), 



X aP (r-r') = -i dtc- rlt ([S a (r,t),S f3 (r',0)}). (A2) 
Jo 



From the conservation of total spin we have x a/3 = 

The spin operators split into forward scattering, Sf(r), 
and backscattering parts, Sg(r). While the forward scat- 
tering contribution to the susceptibility remains regular, 
the backscattering contribution has (at zero tempera- 
ture) a singularity at momenta q = ±2kp. Since this 
singular behavior dominates the physics described in this 
work, we thus keep only the backscattering part. We then 
define 



Xl(r,t) = -i{S?(r,t)S?(0,0)), 



(A3) 



such that Xa(f) is the u> — ► + irj Fourier trans- 
form of X a(r,t) = tf(*)[x>M) - X>(-r,-t)] = 
2#(t)Im[iXoO",i)], with ■& the step function. The opera- 
tors 5" are given by 



(A4) 



tors because they drop out in the averages. We then find 



-j COs(2fc F r) , e i A /g(^ e ( r ,t) + 9.( r ,t)) e -iV2(^ e (0)+6l.(0))v 

(2tto) 2 



(A8) 



Xy(r,t) 



-lCOs(2k F r) / e iv^(^ c (r,t)-g J (r,t)) e -iV2(0e(O)-9 J (O))\ 

(27ra) 2 

(A9) 



= »COs(2fc F r) iV2(^(r,t)+^,(r,i)) e -i^(0 c (O)+^(O))^ 

(27ra) 2 

(A10) 

The determination of these averages is a standard calcu- 
lation in the bosonization technique. At zero tempera- 
ture, wc obtain 63 



-i cos(2kpr) 
(2ira) 2 



r — vpt + iS 



r + vpt — i5 



(All) 



with g x = g y = (K c + K- 1 )/2, g z = (K c + K s )/2, and 
S > a short distance/time cutoff. For g a well below 
1 we can choose S infinitesimal and proceed analytically 
as below. For g a < 1, the singularities in Eq. (|A11[) 
become too pronounced and the cutoff plays an increas- 
ingly important role. We then must include a finite 5 
and proceed numerically with the further calculation be- 
low and the subsequent evaluation of, for instance, the 
cross-over temperature T* (see caption of Fig. [6]). For 
g a < 0.8, however, such a cutoff is not required. For the 
explicit analytic expressions below we keep therefore an 
infinitesimal S. 

The finite temperature expressions are 



with a a the Pauli matrices. Using the bosonization iden- 
tities (|2T)]) and (f2~Tj) wc see that these spin operators can 
be written in the formS S£ = [0% DW + (Og DM /) f ]/2, 
where 



^SDW 



^SDW 



O 



SDW 



p -2ik F r 

2na 

p -2ik F r 

i- e 1 ^ 

2na 

2na 



e iV29 a _ e -iV28 3 



i\p24> s 



(A5) 
(A6) 
(A7) 



We have omitted in the latter expressions the Klein fac- 



-icos(2kpr) I i^ajjivp 
(27ra) 2 \^ sillh ^_£_( r _ UF t 

7ra/ (3vf 



iS 



sin M -^( r + v i pt - iS ) l 



(A12) 



where (3 = 1/kgT. Twice the imaginary part of these 
susceptibilities determines \a i r , t) . We have at zero tem- 
perature 



Xa (r,t) = -#(t)d(vt-\r\) 



sm(-Kg a ) cos(2fc F r) 
27r 2 a 2 



,,2 _ v 2f2 

(A13) 



Sc. 



2G 



and at finite temperature 
Xa (r,t) = -§(t)d(vt-\r\) 



ira/ (3v F 



8m(irg a ) cos(2k F r) 

2~^~ 2 

ira/ f3v F 



sinh ( W^( r + y F t ) 



(A14) 



Further Fourier transforming to (q, uj) space and taking 
the ui — » + irj limit leads to>& 



x a (q) 



sin(7Tff a ) 
Att 2 v f 



E 



K=± 

at zero temperature and 



a(q + n2k F ) 



(A15) 



Xa(q) 



E 



shnV^)^^ ^ s ( f3v F \ 



2-2 9o 



An 2 v F 



-r 2 (i-. 9 a) 



V 2?ra 7 



^(9 + *) 



(A16) 



at finite temperature, where T is Eulcr's Gamma func- 
tion. 



With feedback from the Overhauser field 



We have seen in Sec. IV Al that the feedback from the 
nuclear magnetic field strongly renormalizes the electron 
system and opens a gap B* y for the field 4>+ oc (j) c + 
9 S . Because the fluctuations of this mode are frozen out 
at ksT < B*, y , the response of the electron system to 
external perturbations is modified. In particular, this 
affects the electron spin susceptibility Xa , and here most 
importantly the exponents g 7 which become smaller and 
anisotropic. 

The calculation of Xa in this situation has been carried 
out in Ref. [76| (see also Ref. For completeness we 

outline this calculation here again. 

Let us set f = (r, t), Xa{f) = —i2(na) 2 Xa(r,t). We 
focus first on Xx . Using the relations (pB|) - PT)) and Eq. 
151) we can then write 



Xx (f) = cos(2ft F r)[([e iV33? *+^ , e - 4 ^*+<°)]) 

+ ^ c tV2^<P-(f)-iV2K 7T 4>+{f) ; e iv / 27? 7 0-(O)+iv / 2if 77 + (O)j^ 

(A17) 

with K' = AK C /K S K and K" = (K c - Kj^/K, and 
where we have used the invariance of the Hamiltonian 
P5|) under a simultaneous sign change of all the boson 
fields. 

Since the Hamiltonian (|49)l is quadratic in the boson 
fields all these averages are fully determined by the 2- 
point correlators^ 

7TV — 

{4>*_ (q, uj 0_ (q, uj)) = . — (A18) 

[uj ± ir\Y — vLq 



for the masslcss field and 



(<j)* + (q,uj)(j> + (q,uj)) 



(uj±ir 1 y-v% q *-(B*)* 



(A19) 



for the massive field. The sign of the infinitesimal shift 
Hrj is determined by the time order of the operators. The 
corresponding 9± correlators (important for Xz ) can be 
obtained, for instance, through the equations of motion 
dt<t>±{r, t)=v±V0±(r,t). 

The correct treatment of the singularity in Eq. (|A18[) 
at (u) ± V— q) — > results in the singular power laws (at 
T = 0) of Xa of the LL theory and so leads to a contri- 
bution to the susceptibility entirely equivalent to those 
of Appendix IA 11 

This singularity is absent in Eq. (|A19[) and hence the 
power law singularities are broadened to resonances with 
a width and height determined by B* y . The precise eval- 
uation of the expressions is cumbersome. Approxima- 
tions can be found in Refs. mumm 

For our pur- 
pose, however, this evaluation is not required because the 
physics of the combined electron-nuclear spin system is 
entirely dominated by the singular behavior of the <p- 
correlators at q = 2k F , provided that fc^T < B* y , which 
is a necessary condition anyway that is indeed satisfied 
for the considered systems. 

To see this in detail, let us expand Eq. (|A17p in powers 
of 4>+ ■ The lowest nonzero term in </> + is 

cos(2fc F r) [2Jr([0+(f), 0+(O)]) 

+ 2K"{[e i ^^U+{f) , e- 4VW7 *-(°)V + (0)]) 

(A20) 

The Fourier transform of the first term can be directly 
evaluated with Eq. (|A19[) . In the limit u> — > and 
q — > ±2k F it tends to a constant ~ l/(B* y ) 2 and so con- 
tributes only insignificantly to Xx{q)- The second term 
leads to a sum of expressions of the type 



cos 



{2k F r)e K> <*- ( f )*- <°» (<j> + {f)4>+ (0)) . 



(A21) 



The fields 0_ arc gapless and so lead to power laws of the 
form (|A15|) (replacing q 2 by uj 2 /v 2 _ — q 2 and the exponent 
by 1 — K'). The Fourier transform of the latter expression 
is then a convolution of the form 



dq'duj' 



1 



uj' 2 — v 2 _q a 



l-K' 



[uj-uj') 2 -v 2 + (q ± -q') 2 -{B* xy ) 2 ' 



(A22) 



with q± = q ± 2k F . In the static limit u) — » and at q — > 
±2k F , the uj' integral is dominated by the poles at uj' = 

i^Jv+q' 2 + (B* y ) 2 and the intcgrable weak singularities 

at uj' = ±v_q' . A potential singular behavior must thus 
come from the poles. If we focus on the pole at uj' = 



27 



l v 2 ,2 



(B* ) 2 we obtain 



dq' 



l-if' 



i; 2 g' 2 



(A23) 



The remaining integral leads to an arcsinh, which has 
an ultraviolet divergence that has to be cut off at 1/a. 
More importantly, however, the result has no infrared 
divergence, meaning that this expression remains regular 
at q± — ► 0, with a value determined by B* y . 

From this calculation we conclude that the Fourier 
transform of (|A20[) remains regular at ui — > 0, q± — ► 0. 
Since the Hamiltonian is quadratic in the boson operators 
higher order correlators are products of the latter results 
and so remain fully regular. We have therefore shown 
that the singular behavior of the susceptibility Xx{o) is 
fully controlled by the massless fields 0_ only. This al- 
lows us to entirely neglect the gapped fields and approx- 
imate the susceptibility by 



Xx (f) = cos(2fc F r)( [ e i^2W<f>-(r) 



/2^0_(O)j^ 



(A24) 

which is of precisely the same form as the susceptibility 
of a standard LL as discussed in Appendix IA II The 
difference is the modified exponent K', the velocity u_, 
and an amplitude reduced by 1/2 because the second 
term in Eq. (|A17[) depends on <j) + only and drops out. 
From the results of Appendix IA II and the Hamiltonian 
(|6ip we obtain, for instance, at zero temperature, 



Xx(q) 



1 sm(irg' x 
2 



-r 2 (i-^) 



E 

K=± 



a(q + n2k F ) 



(A25) 

with g' x = K'/2. 

The evaluation of Xy(l) leads to precisely the same 
result, Xv(o) = Xx(q) and in particular g' y =g' x . 

The susceptibility Xz{q) is different in that it involves 
different boson fields. From Eq. (|A7|) we see that Xz 
depends on the combination <fi c + <fi s = ■^=(K c <fi+ — 



K a 

leads to 

Xz(q) = 



+ + ^JK C K S 6-), where we have used Eqs. 
Neglecting the gapped fields </>+ and 9 + then 



. a r 2 l 

4w_7T 2 



g'z) 



E 

K=± 



a(q + n2k F ) 



2-2 s ' 



(A26) 



The prefactor 1/2 of Xx,y is here missing and the expo- 
nent is g' z = {K c /K s + K C K S )/2K. 

The extensions to finite temperatures are obtained in 
precisely the same way as in Appendix IA 11 



APPENDIX B: REAL SPACE FORM OF RKKY 
INTERACTION 

The real-space form J a (r) of the RKKY interac- 
tion can be found by time integrating the susceptibility 



Xa {r, t) , given in Eq. (|A14[) ■ This integration was evalu- 
ated in the limit \r\ -C At in Ref. [24|. Here we determine 
Xa(f) and so J a (r) in the general case. For convenience 
we shall drop the index a = x, y, z in this section. 

The stationary real space form x{ r ) corresponds to the 
lo = Fourier transform of X (r, t), which from Eq. (|A14[) 
is determined by the integral 

sin(7rg) cos(2fcp2;) 



X(r) = — 



27r 2 a 2 
ira/ (3v F 



dt 



\r\/v f 



ira/ f3v F 



(Bl) 

Setting y = tv F /\r\, £ = TT\rU/3v F , and using the relation 
sinh(a+6) sinh(a — b) = sinh (a) — sinh 2 (6) we can rewrite 
this as 

sin(7rg) cos(2k F r) ( ira \ 2g \r\ 
K (3v F J v F 



X(r) 



2ir 2 a 2 



dy I sinh 2 (£) — sinh 2 I 9 . 



(B2) 



A further change of variable y — > z = sinh 2 (£)/ sinh 2 (£?/) 
leads to 

X(r) = - 

r.1 



sin(7n?) cos(2k F r) / 7ra 



2g 



\r\ v F 



2ir 2 a 2 \j3v F J 2£sinh 2ff (£) 

dzz s_1 (l - z)~ 9 [l + z/ sinh 2 (£)]" 1 I 2 , (B3) 

which is the standard form for the Gaussian hypergeo- 
metric function [Ref. [9(| §15.3.1], 

F ( a ' 6;c;z ) = ¥n^r h\ f dzz h -\i- z y- h -\i-wz)- a . 

T(b)T(c - b) J 

(B4) 

Hence, noting that v F /3 = At is the thermal length, 
sin(7r<7) cos(2k F x)\T ( wa/Xx 



^ ^ 4TT 2 a 2 v F \^sinh(7r|r|/AT) / 

x F (1/2, g; 1; - sinh- 2 (^|r|/A T )) . (B5) 

The hypcrgeometric function has the asymptotic expan- 
sion for \w\ -> oo [Ref. H §15.3.7] 

T(b-a)T(c) 



F(a, b; c; w) = 



T(b)T(c-a) 



(-wy^l + Oiw- 1 )] 



T(a-b)T(c) 



(-wy^i + oiw- 1 )]. 



(B6) 



T(a)T(c-b) 

For g > 1/2 we obtain from this expansion the asymp 
totic behavior for |r| -c At [cf. Ref. |24| . 

X(r) 



1 T(g - 1/2) cos( - 2fc ^ ( a 
4irav F T{g)y/ii~ \|r| 
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(B7) 



For g < 1/2 the latter equation must be complemented 
by further corrections from Eq. (|B6[) . The RKKY in- 
teraction is eventually determined by J(r) — A 2 ax(r)/2, 
with r running over the sites of the nuclear spin lattice. 
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APPENDIX C: MAGNON SPECTRUM 

The magnon spectrum about the helical ground state 
of the nuclear spins can be evaluated by mapping the 
helical state back onto a ferromagnct. This allows us to 
use the standard results^ for spin waves in anisotropic 
ferromagnets. We focus on the + hclicity in Eq. (|33|) . 
The case for the opposite helicity is equivalent. 

The mapping is achieved by noting that the ground 
state (|33p can be mapped onto a ferromagnetic align- 
ment along the spin e x direction by defining the local 
transformation I4 = with 



cos(2fci??"i) sin(2feF7"i) 
R%= I — sm(2kpri) cos(2k F ri) 
1, 



(CI) 



In the classical ground state we then have I; = 
(/JVj.,0,0). The RKKY Hamiltonian (13|) becomes in 
terms of these new spin vectors 



H 



off 



a/3 



zj.a/3 J ~ 







with 3?f = ^RTRfJl- Hei "e the RKKY couplings 
are generally assumed to be anisotropic as J?. = Jf^ ^ 
, which covers both the isotropic case without the feed- 
back, and the anisotropic case with the feedback on the 
electron system. In Fourier space the latter Hamiltonian 
becomes 



1 

N 



E 

9^0, a/3 



jx 



Nl 



R 



a/3 
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Nl 
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a/3 



q ry 
N 2 ^0 
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(C3) 



(C4) 



The magnon description is obtained by replacing the 
spin operators I" with the Holstein-PrimakolSi bo- 
son operators a q . This leads to I q = (IN±)NS q fl — 
jrZ P 4+ g <h>, h = (IN ± /2)^(al q + a q ), and JJ = 
— i(IN_i/2) 1 / 2 (a_ q — a q ). Neglecting terms involving 
more than two a q operators results in a Hamiltonian that 
is equivalent to (|C3j) up to corrections of order 1/IN±. 
Since Nj_ 3> 1 this approximation is very accurate. We 
then obtain 



H — Eq 



-Y^{al,a- q ) 

q>0 





h q 2) ) 








fe) 



with Eq 
(-2^ 



JX 

J 2k F -q 



(N ± I)(N ± I - 
J z q )l N - 



m(J2k F /Nl), h, 
(J, 



(1) 



and h, 



(2) 



2k- T 



J q )/Nj_. Finally diagonalizing this Hamiltonian leads 
to two magnon bands with dispersions 



(i) 
1 

(2) 



= 2i{j% kF _ q - J% kF )/N ± , 



2I(J*-J% kF )/N 



(C6) 
(C7) 



The first magnon band is gapless and dominates the fluc- 
tuations about the nuclear ground state. The second 
magnon band is gapped at q = with a gap of the 
size ~ I Jjj, I ■ Since quite generally through the feedback 

(2) 

|J| feF | > IJfkjJ we have, however, u q ' < at q ~ 2fc^. 
This normally leads to an instability of the ground state 
and so normally means that the assumption of the or- 
dered ground state is not valid. However, as shown in Sec. 
IVI C[ in the present case such a destabilization would de- 
stroy the feedback that stabilizes the order and so causes 
the anisotropy J q ^ J q . This destruction would have a 

very high cost in ground state energy. Therefore the oc- 

(2) 

cupation of the u q < modes by a macroscopic magnon 

(2) 

(C2) number is not possible. The remaining u> q > modes 
can then be captured by the same treatment as the ujq 1 ^ 
modes. Yet as they involve large momenta ~ 2kp , they 
are neglected in the present treatment. 



APPENDIX D: MAGNETIZATION FOR 
SYSTEMS WITH SMALL LEVEL SPACING 

For systems with a large enough length L such that the 
consistency relation ^> fesTg (or T*) put forward in 
Sec. II V CI no longer holds, we need to reexamine carefully 
the derivation of Eqs. ([55)1 and (|70|) for the magnetiza- 
tion vnik F ■ This situation is indeed met for typical GaAs 
quantum wires (see Table H]) where A^ ~ ksT^. In this 
case the singularity in the magnon occupation number in 
Eq. ([55)) may be controlled by the cutoff at q = ir/L, 
i.e. by Ai instead of the L-independent scale /cbTq 
(or ksT*). Here we show, however, that this singular- 
ity dominates the magnetization and hence modifies the 
cross-over temperature only for L that lie many orders 
of magnitude above any realistic length. The result (|70p . 
therefore, remains valid and robust also for A^ > /cbTq*, 
and in particular also for A^ ~ k B T* . 

To see this we focus on the q — > part only in the 
magnon occupation number in Eq. (|35[) and neglect the 
contributions u q ~ 2/| J 2 fc F |/iVj_ that have led to Eqs. 
(f3"8")) and (|TD|) before. For simplicity we discuss only the 
case without feedback, and note that the conclusion is the 
same for the case with the feedback. For q < tt/\t we 
have uj q < 2I\J2k F \/N±_ and the dispersion u q can be well 
approximated by a parabola, ujq/ksT = Dq 2 , with D = 
{kBT^dPujq/dq 2 ^ = 2ICv 2 F (k B T) 2 9- 5 /N ±: and C = 
(C5) C(g,v F )T 4 (g/2)[7T 2 - 2sin 2 (7r 5 /2)^ 1 ( 5 /2)]/16 7 r 4 . Here 
ipi is the polygamma function. 
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We then ask under which conditions the following ap- 
proximation holds: 



i 



2a 



IN ± N ^ e^ q /k B T _ 1 

q>0 



IN, N 



7T/L 



dq 



1 



2?r e D <i 2 - 1 ' 



(Dl) 

At T < T * (T*) the integrand at q = n/ \ T is small and we 
can push the upper integral boundary to infinity. Setting 
then x = \J~Dq and X = \/Dit/L = Lu^/L/ksT we have 
to evaluate 



Q(x) 



dx 



x 

dx 
x x 
1.294 



1 Jo 

fX 

dx 
'o 

1 X 

x + Y 



dx 



1 



1 



1 



1 



c x — 1 X 1 
2\ 



0{X 2 ). 



(D2) 



At the considered temperatures T < T *, X is a small 
number, and we can keep the 1/X contribution only. 
This leads to 



where 



k B T* L = [CI 2 TT 2 /La\ 



(D5) 



is an L-dependent characteristic temperature. This tem- 
perature defines the cross-over temperature if T£ < T * . 
Putting in numbers corresponding to GaAs quantum 
wires (sec Table Q} we see, however, that T£ > T *. For 
T < Tq, therefore, the T£ dependent contribution to 
Tn2k F is negligible. Only at entirely unrealistic lengths 
L ~ 1 m wc obtain T£ ~ T Q *. For any realistic sam- 
ple, therefore, T£ and so L do not affect magnetization 
and cross-over temperature even for Al < ksT^ and Eq. 
(|38|) remains valid. As stated above the same holds for 
the case with feedback where Tq is replaced by T* and 
the magnetization is given by Eq. ([701) . 



m 2 k f 



1 - 



2a 
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ILN ± VD(tt/X t ) y/u*, L /k B T 



2ak B T 
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